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ABSTRACT 

The  decision  problem  of  the  allocation  of  available  weigt 
space  in  an  air-dropped  anti-submarine  torpedo  to  fuel^  explo&i 
to  attack  speed  capability  above  some  preset  minimum  is  studied  for 
case  in  which  alternatives  under  consideration  do  not  differ  with  respect 
to  their  effect  on  weapon  reliability  and  to  their  effect  on  enemy 
countermeasures  capability „  Analytical  models  are  developed  which 
relate  (1)  probability  of  acquisition  to  endurance  for  the  case  of 
circling  search  at  constant  depth,,  (2)  probability  of  hit  given  acqui- 
sition to  endurance  and  attack  speed,  (3)  probability  of  hit  on  a 
close- in  attack  run  to  attack  speed,  and  (k)   probability  of  kill  damage 
given  a  hit  to  amount  of  explosive  used. 
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TABLE  OF  SYMBOLS  AND  ABBREVIATIONS 
A.  =  area  in  the  j  th  plane  which  is  scanned  during  one  ping  cycle 


J 


,2 


(yd  );  j  =  J',  ••,  J 


?  i 


A.,(B)  =  echo  excess  in  plane  j  at  a  range  of  100k  yards  which  is 
J*1 

required  to  acquire  a  target  at  an  azimuth  angle  (  p")  to 

the  direction  of  the  transmitted  ping  (db)j  J  ~  ~3S    °  °  °  $   J| 

AMV,   =  effective  volume  in  level  L,  above  weapon  search  depth  but 

below  the  critical  depth  missed  inside  the  insonified  volume 

on  ping  zero  on  the  u  th  search  turn  (yd  )j 

_  — • 

d  =  D  +  1,  ••  •,  D  ;  u  =  1,  •  ••,  U 

ASV    =  effective  volume  in  level  L  above  weapon  search  depth 
scanned  on  ping  zero  on  the  u  th  search  turn  (yd  )j 
d  ■  1,  •••,  D  ;   u  ■  2g  •••,  0 

B.,     =  echo  excess  in  plane  j  at  a  range  of  100k  yards  available 
J* 

for  acquisition  on  the  azimuth  reference  axis,  (db) 

j  =  -  J,  —,  J;  k  ■  1,  •••,  Kj 

BMV,   =  effective  volume  in  level  L,  below  weapon  search  depth  but 
above  the  critical  depth  inside  the  insonified  volume  on 
ping  zero  on  the  u  th  search  turn  (yd  ) 
d  =  D  +  1,  °",    D1  |  u  a  1,  •••,  U 

BSV,   =  effective  volume  in  level  L,  below  weapon  search  depth 
scanned  on  ping  zero  on  the  u  th  search  turn  (yd  )j 
d  =  1,  •••,  D*  j  u  =  2,  ...,   U 

BW     =  length  of  the  beam  dimension  of  the  target  (yd) 

b.     =  normalized  pattern  function  in  the  vertical  plane  for 

transmitting  pattern  (a  function  of  vertical  off=axis  angle) 

b      =  normalized  pattern  function  in  the  vertical  plane  for 

receiving  pattern  (a  function  of  vertical  off-axis  angle) 
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b_     =  normalized  pattern  function  in  the  azimuth  plane  for  trans<= 
mitting  pattern  (a  function  of  azimuth  off-axis  angle) 

b.     =  normalized  pattern  function  in  the  azimuth  plane  for  receiving 
pattern  (a  function  of  azimuth  off-axis  angle) 

CDA    =  critical  depth  above  search  depth  above  which  a  target 

trapped  inside  the  insonified  volume  may  escape  and  below 
which  there  is  no  escape  (ft) 

CDB    =  critical  depth  below  search  depth  below  which  a  target 

trapped  inside  the  insonified  volume  may  escape  and  above 
which  there  is  no  escape  (ft) 

CH     =  distance  between  the  weapon  and  the  intersection  of  a  ping 
bisector  and  the  target  cylinder  radius  at  the  time  of  ping 
transmission  (yd) 


CR.,    =  distance  from  the  center  of  the  target  cylinder  to  the 
projection  of  the  end  point  of  radial  (R..  5  P  ..)  in  the 
horizontal  plane  (yd);  j  -   J°  ,  ■••,  j";  k  s  1,  »•<>,  K, 

c      =  velocity  of  sound  in  the  medium  (yd/sec) 

D,  D   =  value  of  the  index  d,  respectively  below  and  above  search 
depth, when  the  boundary  nearest  search  depth  of  the  level 

L,  is  the  critical  depth  (an  integer) 

i  _i 

D  ,  D  =  upper  limits  on  the  numbers  of  levels  L,  respectively  below 

and  above  search  depth  (an  integer) 

D      =  minimum  target  depth  to  be  considered  (ft) 

D_-.    =  maximum  target  depth  (ft) 
in 

D      =  weapon  search  depth  (ft) 
w  & 

d      =  index  of  the  levels  L,  above  and  below  search  depth  (an 

-•  ° 

integer);  d  =  1,  •••,  D  and  d  -  1,  •••,  D  respectively 

F(w)   =  In  P,  (w),  (a  negative  real  number) 
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f'(w)  ■  F(w)  -  Xk-Cw^)  +  g  (w  )  +  g  (w  )  |,  (a  negative  real  number) 

'  d(x)  s  cumulative  normal  probability  distribution  with  mean  zero 
and  standard  deviation  one  evaluated  at  the  argument 
x  (a  decimal) 

f  (x)  =  probability  density  function  of  the  random  variable  hit  depth 
HD 

evaluated  at  the  argument  x  (a  positive  real  number) 
f  (w.  )  =  In  P  (w1  ),  (a  negative  real  number) 

f2^W2jk'  W2ik^  =  ln  PHjk^W2Jk'  W2jk^  ^a  neSative  real  number)^ 

i-1    co.    T  «  W  —  1    000   If 

f  (w  )   =  In  P  (w_)  (a  negative  real  number) 
5     j         E  j 

G,     =  function  relating  increased  weapon  weight  to  increased 
weapon  speed  above  S  (a  function  of  w   ) 

G_     =  function  relating  weapon  speed  and  depth  to  fuel  expenditure 
rate  (a  function  of  speed  and  depth) 

g.(w.)  =  function  relating  component  weight  to  volume  (a  function  of 
w^  i  =  1,  2,  2',  2",  3) 

g.     =  P  ("acquisition  during  the  i  th  ping  cycle  J  ,  (a  decimal)^ 
i  =  1,  •••,  I 

g.     =  P  [acquisition  during  the  i  th  ping  cycle  of  the  u  th  search 
turn  J  ,  (a  decimal)^  i  =  l9    •  ••,  M^  u  =  lp  »•»,  0 

H.     =  altitude  of  an  isosceles  triangle  of  area  A.  (yd)j  J  ■  Jc 9 
"-,    J" 

H,     =  one  way  transmission  loss  at  a  range  of  100k  yards  (db)* 
k  =  1,  ••»,  Kj 

HD     =  depth  at  which  the  weapon  hits  the  target^,  a  random  variable  (ft 

HD     =  average  hit  depth  (ft) 

HV  (w_)  =  hull  vulnerability  for  kill  damage  of  the  p  th  target  surface 
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area  element  at  the  midpoint  depth  of  the  q  th  depth  sub= 
interval  given  an  explosive  weight  w  (one  or  zero); 
p  =  1,  •••,  P;  q  =  1,  ...,  Q 

h.(w)    =  recurrence  relation  for  the  allocation  model  evaluated  at 
the  argument  W  (a  negative  real  number);  I  m  X9  29   3 

h       =  P  [hit  on  the  n  th  attack  (or  reattack)  \   s    (a  decimal); 
n  =  1,  ••',  N 

I       =  upper  limit  on  the  number  of  ping  cycles  (an  integer) 

i        =  index  of  the  total  number  of  ping  cycles  (i  -   15  •••,  l) 
of  the  number  of  a  ping  in  the  first  search  turn 
(i  =  1,  •••,  M),  (an  integer) 

J       =  a  limiting  angle  in  the  vertical  plane  for  off-axis  signal 
reception  measured  from  the  transducer  axis  (deg) 

j' ,  J1 '   =  smallest  and  largest  angles  respectively  in  the  vertical 
plane  at  which  acquisition  may  occur  (deg) 

J,       =  the  number  of  combinations  of  w  .  and  w  ,  considered  which 
sum  to  w  ,  (an  integer);  k  ■  1,  <"•,  K 

j        =  index  of  the  planes  perpendicular  to  the  vertical  plane 

containing  the  transducer  which  are  available  for  acquisition 
(j  =-J>  "*t   «0  an(*  *-n  which  acquisition  may  occur 
(j  =  j' ,  *",  j"),  (an  integer);  also  index  of  the  combina- 
tions of  w„..  and  w  '   considered  which  sum  to  w_,  (an  integer) 
2jk      2jk  2k  *  ' 

j  =  1,  •  •  o,  Jfc  k  =  1,  •••,  K 

K       =  number  of  values  of  w  ,  to  be  considered  (an  integer) 

K.       aa  the  number  of  radials  in  plane  j  used  to  outline  the  insonified 
area  (an  integer);  j  =  J*  >  *",  J11 

K1       =  the  number  of  ranges  in  the  plane  j  which  have  a  non-negative 
available  echo  excess  (an  integer)  j  =  -J^  •••,  J 


IADV 


du 


LBDV 


du 


M 
MAI 


du 


MBDV 


du 


M(BG). 


MRH 


m 


=  index  on  the  weapon  payload  weight  used  in  event  2  (k  a  19 
on  the  ranges  in  the  plane  j  under  consideration  (k  a  ls    ° ■ ° s 
and  on  the  radials  outlining  insonified  area  in  plane 
j(k  =1,  •••,  K.),  (an  integer) 

=  effective  volume  element  in  level  L,  at  or  above  the  critical 

d 

depth  above  search  depth  on  ping  zero  on  the  uth  search  turn 


K 


(yd3);  d  =  1,  -..,  D;  u  =  2, 


U 


=  effective  volume  element  in  level  L,  at  or  below  the  critical 

d 

depth  below  search  depth  on  ping  zero  on  the  u  th  search  turn 
(yd3);  d  =  1,  —,  D;  u  =  2,  •••,  U 

=  identification  of  the  d  th  level  (or  layer)  above  ano  below 
search  depth  (d  =  0,  •••,  D  and  d  *  0,  •••,  D  respectively ))| 
also  the  length  of  the  d  th  level  above  and  below  search 
depth  in  the  vertical  cross=section  of  the  first  turn  insonified 
volume  of  ping  zero  (d  =  1,  •••,  D  and  d  s  1,  °  °  °  $   D 
respectively),  (yd) 

=  number  of  ping  cycles  in  a  search  turn  (an  integer) 

=  effective  volume  element  in  level  L,  at  or  above  the  critical 

d 

depth  above  search  depth  on  ping  zero  on  the  u  th  search 
turn  (yd3);  d  =  1,  •••,  Dj  u  =  2,  »••,  U 

=  effective  volume  element  in  level  L,  at  or  below  the  critical 

d 

depth  below  search  depth  on  ping  zero  on  the  u  th  search 


turn  (yd  );  d  =  1,  ■••,  D;  u  = 


D 


=  background  masking  level  at  the  weapon  transducer  (db) 
=  distance  traveled  by  the  weapon  from  time  t  to  time 

'cPA  <*> 
=  index  indicating  the  leftmost  radial  (R.lk^  P-iu)  or  the 

rightmost  radial  (R  .  (   Piou)  ^or  acquisition  at  a  range  of 

100k  in  plane  j  (one,  two) 
=  total  number  of  attacks  (or  reattacks)  completed  on  the 

target  (an  integer) 
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n       =  index  for  the  attacks  (or  reattacks)  on  the  target;  (an 
integer)  ;  n  =  1,  •••,  N 

P        =  the  number  of  target  surface  area  elements  (an  integer) 

P  [x]     =  probability  that  event  x  occurs  (a  decimal) 

P       =  P  [target  acquisition  (event  1)JP  (a  decimal) 

P       =  P  [kill  damage  is  inflicted  on  the  target  by  explosion  of 

the  weapon  warhead  (event  3)  given  that  event  2  has  occurredj 
(a  decimal) 

P       =  P  [hit  occurs  on  the  p  th  target  surface  area  element]  s 
(a  decimal);  p  =  1,  ««.,  P 

P_       =  P  fhit  occurs  at  the  midpoint  depth  of  the  q  th  depth 
Eq  fc 

subintervalj  ,   (a  decimal);  q  -  1^  °>°9   Q 

P       =  P  [target  is  hit  (event  2)  given  that  event  1  has  occurredj 9 
(a  decimal) 

P  ,     =  P  for  the  j  th  combination  of  w'  ,  and  w°°   making  up  the 

HjK  H  tJK       cJK 

weight  w  ,  for  event  2;  j  =  1^,  °  °  °  9    J,|  k  =  1,  •••,  K 

P.       =  P  [target  kill]  ,  (a  decimal) 

PQ.      =  P  [target  is  in  the  i  th  volume  slice  of  the  target  volume 
cylinder  at  the  commencement  of  searchjp  (a  decimal); 
i  =  1,  ■••,  M  -  1 

PRR      =  weapon  pulse  repetition  rate  for  search  (pings/sec) 

p       =  index  of  the  target  surface  area  elements  (an  integer); 
p  =  1,  ...,  P 

Q       =  number  of  hit  depth  subintervals  (an  integer) 

q       =  index  of  the  hit  depth  subintervals  (an  integer); 

q  m    1,  ...,  Q 

R       =  number  of  values  of  w.   to  be  considered  (an  integer) 
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R 


RA 

RAD\ 


RBDV 


R, 


R. 


=  horizontal  distance  from  the  weapon  to  the  center  of  the 
target  volume  cylinder  (yd) 

=  P  target  kill  ,  (a  decimal) 

=  acquisition  range  (yd) 

=  effective  volume  element  in  level  L,  at  or  above  the 
du  d 

critical  depth  above  search  depth  on  ping  zero  on  the 
u  th  search  turn  (yd  );  d  =  1,  ° ° •,  Dj  u  -  2$    ,,oj  0 

,    =  effective  volume  element  in  level  L,  at  or  below  the 
du  d 

critical  depth  below  search  depth  on  ping  zero  on  the 

o 

u  th  search  turn  (yd  );  d  .2  1,  °  •  ° $   D;  u  as  28    •••,  U 

=  radius  of  the  target  volume  cylinder  at  the  time  the 
weapon  splashes  into  the  water ,  (yd) 

=  radius  of  the  target  volume  cylinder  at  the  start  of  the 


i  th  ping  (yd);  1  =  1, 


M 


Ijk 


maximum  acquisition  range  at  angle  P   in  plane  j  on 
i  th  ping  (yd);  i  =  1, 
•",  K. 


M  -  lj  j  =  J"«  -"J" 


kjk 


r' 


jk 


Ljmk 


RH 

r 
r 


=  length  of  the  k  the  acquisition  radial  in  plane  j  (yd 


J  =  A 


,  j";  k  as  lt 


K 


=  length  of  the  projection  in  the  horizontal  plane  of  the 
k  th  acquisition  radial  in  plane  j  (yd)|  j  «  J° ,  ■•», 

j"  ;  k  =  1,  ....  K. 

J 

=  length  of  the  acquisition  radial  (R,  ,  j  pe  ,  )  in  plane  j 
(yd);  j  =  J'  $    •",  J";  ni^  1,  2;  ks  ls  « ••,  K° 

=  range  from  the  weapon  to  the  target  center  at  time  t  (yd) 

=  index  of  the  weapon  payload  weight  used  in  event  1  (an  integer 

=  coordinate  of  the  intersection  of  the  weapon  track  extension 
with  the  vertical  plane  (a)  through  the  submarine  target 
center  perpendicular  to  any  plane  containing  both  target 
and  weapon  at  time  t  (yd) 
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r1       =  coordinate  of  the  intersection  of  the  weapon  track  with 

the  vertical  plane  (A*  )  parallel  to  plane  A  containing  bon 
the  submarine  target  center  and  weapon  at  time  t    (yd) 

S       =  number  of  values  of  w_   to  be  considered  (an  integer) 

3s 


S        =  preset  minimum  value  of  attack  speed  (yd/sec) 
o 

S  =  target  cruise  speed  (yd/sec) 

S_w  =  target  maximum  speed  (yd/sec) 

S  =  target  speed  during  weapon  search  (yd/sec) 

S  .  =  weapon  attack  speed  (yd/sec) 

S  =  weapon  search  speed  (yd/sec) 

SL.,     =  signal  level  of  the  target  echo  at  range  100k  in  p 
(db);  j  =  -J,  -..,  J;  k  s  1,  ...,  K^ 

SV.      =  effective  scanned  volume  on  the  i  th  ping  cycle  of  the 

u  th  search  turn  (yd-*)j  i  -  1^,  '••,  M  =  lj  u  -  1,  °  °  ° $   U 

s        =  index  of  the  weapon  pay load  weight  used  in  event  3  (a«  integer 
s  =  1,  •  •  •,  S 

T       =  time  used  by  weapon  at  attack  speed  during  reattack  (sec) 

T       =  time  used  by  weapon  at  search  speed  during  reattack  (sec) 
s 

TES      =  distance  traveled  by  target  during  one  complete  search 

turn  (yd) 

TJ       =  index  value  of  the  reference  plane  j  used  to  define  level 
L  in  the  calculation  of  the  special  volume  elements  (an 
integer) 

TS       =  target  strength  (db) 

TV.      =  volume  of  the  i  th  volume  slice  of  the  target  cylinder  at 

i     3\ 
commencement  of  weapon  ping   (yd  )°    i  -  1,  °  °  ° »   M  •=  1 
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t       =  the  time  weapon  sensors  lose  contact  with  the  target  on  an 
attack  or  reattack  (zero  reference  time) 

t       =  time  from  t  at  which  the  weapon  completes  the  reattack  turn 
and  reacquires  the  target  (sec) 

t       =  time  from  t  of  closestpoint  of  approach  of  weapon  and  targe* 

CrA  O 

(sec) 

t„      =  dead  time  from  t  at  which  the  weapon  commences  the  reattack 
DT  o  r 

turn  (sec) 
U       =  maximum  number  of  search  turns  (an  integer) 
u       =  index  of  the  weapon  search  turns  (an  integer)^  u  ~   1,  •••,  U 

V  =  available  payload  volume  (ft  ) 

V  =  j  th  volume  element  of  the  insonified  volume  on  ping  zero 

(yd3);  j  =  j',  ...,  JM 

v.       =  fuel  volume  for  acquisition,  event  1  (ft  ) 

v      =  payload  volume  for  hit,  event  2  (ft  ) 

v'     =  fuel  volume  for  hit,  event  2  (ft  ) 

v' '      =  volume  for  modifications  to  the  propulsion  system  to 

provide  a  maximum  (attack)  speed  capability  in  excess  of 

S  (ft3) 

o     ' 

v0.,     =  j  th  value  of  payload  volume  for  event  2  when  payload  weight 
2  jk 

for  event  2  is  w2k;  j  =  1,  •••,  Jfc;  k-  1,  •••,  K 

v'  ,     =  j  th  value  of  fuel  volume  for  hit,  event  2,  when  payload 
2jk  _ 

weight  for  event  2  is  w2k  (ft  );    j  =  1,  ooo^  Jfc;  k  =  1,  >  ° « » 

v' '      =  j  th  value  of  volume  for  attack  speed  capability  when  p 
2jk  _ 

weight  for  event  2  is  w  ,  (ft  )j  j  -  1,  «■>»,  J,;  k  ~  1^ 


v-       =  explosive  volume  (ft  ) 

W       =  available  payload  weight  (lb) 
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t   _l 
W  ,  W    =  lower  and  upper  limits  respectively  on  the  allowable  range 

of  fuel  weight  (lb) 

W  ,  W   =  lower  and  upper  limits  respectively  on  the  allowable  range 

of  weight  for  modifications  to  the  propulsion  system  to 

provide  a  maximum  speed  capability  in  excess  of  S  (lb) 
mi   _i  1 1 
W   ,  W    =  lower  and  upper  limits  respectively  on  the  allowable  range 

of  explosive  weight  (lb) 

WEL,  =  weapon  enabling  level  at  a  range  100k  yards  (db) 

WS  =  weapon  on-axis  source  level  (db  at  1  yard  from  transducer) 

w.  =  fuel  weight  for  acquisition,  event  1  (lb) 

w.  =  r  th  value  of  acquisition  fuel  weight  (lb);  r  -  1,  »••,  R 

w  =  payload  weight  for  event  2  (lb) 

w  ,  =  k  th  value  of  payload  weight  for  event  2  (lb)j  k  ~  1^,  ooo^  K 

w  =  fuel  weight  for  hit,  event  2  (lb) 

w       =  weight  for  modification  to  the  propulsion  system  to  provide 

a  maximum  (attack)  speed  capability  in  excess  of  S  (lb) 

o 

w  ,     =  j  th  value  of  fuel  weight  for  hit,  event  2.s   when  payload 

weight  for  event  2  is  v^   (lb);  j  ^  lp  »•,  Jfc;  k^  1^  .«.,  1 

w  ,     =  j  th  value  of  weight  for  attack  speed  capability  when 


2jk 


payload  weight  for  event  2  is  w2k  (lb)j  j  s  1,  ...,  j^j 
k  =  1,  •••,  K 


w_       =  explosive  weight  (lb) 

w       =  s  th  value  of  explosive  weight  (lb);  s  s  1,  •»•,  S 

w_       =  average  weapon  fuel  expenditure  rate  in  dive  to  search 
depth  (lb/sec) 

w.       =  weapon  fuel  expenditure  rate  at  search  depth  and  speed  f 

w       =  weapon  fuel  expenditure  rate  at  attack  speed  and  at  the 
average  depth  for  the  attack  (lb/sec) 

xvi 


w       =  weapon  fuel  expenditure  rate  at  attack  speed  and  at  the  mid- 
point of  the  q  th  depth  interval  ( lb/sec )j  q  -  1^  .».,  Q 

w_       =  weapon  fuel  expenditure  rate  at  search  speed  and  at  the  mid- 
Stq 


^S 


o  o 
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point  of  the  q  th  depth  interval  (lb/sec)^  q  s  1 
wn_      =  fuel  used  in  dive  to  search  depth  (lb) 

wl       =  fuel  used  in  the  n  th  attack  or  reattack  (lb) 

2n  %      * 

«  ,«'     =  limiting  angles  for  off=axis  reception  in  the  azimuth  plane 
to  left  and  right  respectively  of  the  weapon  tran 
longitudinal  axis  at  time  of  reception  (deg) 

p'       =  angle  in  plane  zero  the  right  hand  edge  of  a  ping 

the  reference  axis  after  the  overlap  correction  (deg) 

p\,      =  azimuth  angle  of  the  k  th  radial  in  plane  j  (deg)j 

\    •      T*    ...    I  '  •  k  -  1    „  „  „   K 

8       =  azimuth  angle  of  the  acquisition  radial  (R  ,  9    p\  .  )  in 
plane  j  (deg)j  j  =-J,  » .  °  9   Jj  m  ■  1,  2j  k  =  1,  . ..,  K 

Aw  =  increment  of  weight  of  any  of  the  weight  variables  (lb) 

b  =  angle  through  which  the  weapon  turns  between  pings  (deg 

6  '  =  angle  assigned  to  the  M  th  ping  (deg) 

X  =  Lagrange  multiplier  (non- negative  real  number) 

u     =  weapon  angle  of  attack  (pitch)^  (deg) 
Of     O  *  =  standard  deviations  of  the  normally  distributed  random 

variables  r  and  r*  respectively 


0        =  angle  of  dive  to  search  depth  (deg) 

lid 

0        =  angle  through  which  the  weapon  turns  while  a  ping  travel 
to  a  range  of  100k  yards  and  returns  (deg) 

01.       =  weapon  reattack  turn  rate  (deg/sec) 

=  weapon  search  turn  rate  (deg/sec) 
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CHAPTER  I 
INTRODUCTION 
The  Problem 

An  air-dropped  anti-submarine  torpedo  may  be  characterized  by  a 
long  list  of  parameters .   Specification  of  parameter  values  in  a 
specific  weapon  is  made  such  that  the  weapon  may  achieve  a  maximum 
of  effectiveness  in  an  actual  tactical  situation  in  wartime „ 
state  of  the  art  in  weapon  design  will  determine  the  values  of 
of  the  weapon  parameters.  For  example^,  the  transducer  parameters 
may  be  prespecified  by  the  status  of  some  continuing  torpedo  sonar 
development  program.  Cost  considerations  may  also  be  of  great 
importance  in  the  evaluation  of  alternative  proposals  for  the 
design  parameter  values.  For  example^  a  choice  between  two  propul- 
sion systems  may  hinge  on  a  difference  in  development  costs  necessary 
to  produce  a  reliable  weapon.  However^,  for  torpedoes  and  more 
especially  for  air-dropped  torpedoes^  the  restrictions  of  weight 
and  of  physical  dimensions  (or  volume)  resulting  from  the  varied 
mix  of  capabilities  of  the  prospective  weapon  delivery  systems  are 
also  very  significant  in  the  weapon  design „ 

Consider  the  restraints  of  weight  and  volume.  One  of  the  design 
parameters  which  has  a  great  effect  on  weapon  weight  and  volume  is 
the  maximum  operating  depth.  An  increase  of  a  few  hundred  feet  in 
the  operating  depth  requirement  will  require  a  more  substantial 
structure.   This  in  turn  may  result  in  a  decreased  weight  and  volume 
to  be  made  available  for  other  weapon  components „  However^  intelli= 
gence  reports  or  estimates  of  the  depth  capabilities  of  prospective 


submarine  targets  in  the  time  period  for  which  the  weapon  is  to  be 
designed  may  prompt  the  decision  maker  to  fix  the  design  operating 
depth  as  an  input  to  the  design  problem*   Similarly  the  advance  speci- 
fication of  a  minimum  attack  speed  capability  may  be  prompted  by 
intelligence  estimates  and  will  certainly  act  to  place  an  upper  limit 
on  weight  and  volume  available  for  other  weapon  components .   Suppose 
that  a  study  of  specifications  from  higher  authority^,  engineering 
capabilities  to  produce  the  various  weapon  components^,  cost  considera= 
tions  or  some  appropriate  optimization  procedures  have  resulted  in 
specification  of  values  for  all  weapon  design  parameters  except 

1.  Amount  of  fuel 

2.  Amount  of  explosive 

3.  Maximum  speed  (above  some  preset  minimum  value  )„ 

Further,  suppose  the  restraints  on  weight  and  volume  have  been  deter° 
mined  such  that  if  the  amount  of  fuel  and  explosive  and  the  maximum 
speed  capability  above  some  specified  minimum  are  considered  as  a 
weapon  payload,  then  the  available  payload  weight  and  volume  are 
provided  to  the  decision  maker.  The  decision  problem  of  allocation 
of  available  payload  weight  and  volume  to  amount  of  fuel;,  amount  of 
explosive,  and  to  maximum  speed  capability  above  some  specified 
minimum  is  the  problem  considered  in  this  paper 8 


The  Approach 

An  initial  approach  to  the  allocation  problem  is  to  set  certain 
minimum  limits  on  the  parameter  values „  The  minimum  value  of  maximum 
(attack)  speed  has  been  specified ,   There  must  be  sufficient  fuel  to 
provide  some  minimum  probability  of  target  acquisition.,  (in  the  case 
of  continuous  circular  search  at  a  constant  preset  depth  as  is  con- 
sidered in  this  paper  a  minimum  probability  of  acquisition  might 
result  from  a  dive  to  search  depth  and  completion  of  one  complete 
search  turn.)  Then  there  must  be  sufficient  fuel  for  at  least  one 
attack  under  some  assumed  conditions  of  acquisition  range  and  target 
course  and  speed.  Finally  there  must  be  at  least  enough  explosive 
to  result  in  kill  damage  if  a  hit  occurs  at  the  most  favorable  depth 
and  point  of  contact  on  the  target  surface.  Having  determined  these 
allocation  minimums  there  remains  the  problem  of  allocating  the 
balance  of  available  pay load  weight  and  volume „ 

A  second  step  may  determine  a  maximum  value  of  attack  speed „ 
This  maximum  may  be  determined  from  a  study  including  the  probable 
operating  depth,  the  characteristics  of  the  guidance  system^,  and  the 
generated  self-noise  level  and  its  effect  on  the  ability  of  the  sonar 
to  maintain  contact  with  the  target  on  the  attack  run. 

What  is  the  next  step?  How  does  the  decision  maker  decide  between 
enough  additional  fuel  for  one  reattack  and  enough  additional  explosive 
to  produce  kill  damage  over  a  given  additional  element  of  target 
surface?  What  measure  of  effectiveness  is  used? 

This  paper  proposes  a  procedure  which  may  provide  the  decision 
maker  with  a  quantitative  analysis  or  comparison  of  the  tradeoffs 
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involved  in  alternatives  of  allocation  among  finely  explosive,, 
additional  attack  speed  capability.  The  tradeoffs  are  measured  in 
terms  of  probability  of  kill  as  a  measure  of  effectiveness 
ity  of  kill  is  dependent  upon  weapon  reliability  and  upon 
capability  of  the  enemy  to  employ  countermeasures  (in  addition  to 
maneuvering).   It  is  assumed  in  this  paper  that  weapon  reliability 
and  enemy  capability  to  employ  countermeasures  (in  addition  to  maneu<= 
vering)  are  both  independent  of  the  amount  of  fuel  and  explosive  and 
of  the  weapon  attack  speed  within  the  range  of  values  between  the 
preset  minimum  and  some  determinable  maximum.  This  enables  removal 
of  reliability  and  enemy  countermeasures  capability  (less  maneuvering 
from  the  problem  and  probability  of  kill  may  then  be  factored  into 
three  independent  probabilitiesj  those  of  probability  of  acquisition^ 
probability  of  hit  given  acquisition,,  and  probability  of  kill  given 
hit.  An  optimization  is  made  of  probability  of  kill  as  a  function 
of  the  fuel,  explosive,  and  additional  attack  speed  capability  with 
weight  and  volume  restraints.  This  does  not  provide  a  number  for 
use  as  an  operational  value  of  kill  probability.  However  the 
dynamic  programming  methods  used  provide  a  solution  which  specifies 
the  optimal  allocation  of  weight  and  volume. 

It  is  emphasized  that  a  quantitative  comparison  of  the  tradeoffs 
involved  in  the  allocation  process  may  be  obtained  only  by  some 
common  analysis  of  the  three  study  variables  as  they  act  upon  the 
events  of  acquisition,  hit,  and  kill. 


The  dynamic  programming  optimization  model  developed  generates 
an  input  requirement  for  (1)  probability  of  acquisition  as  a  function 
of  available  fuel,  (2)  probability  of  hit  as  a  function  of  available 
fuel  and  attack  speed  and  (3)  probability  of  kill  damage  as  a  function 
of  explosive  weight.  These  requirements  may  be  satisfied  in  a  number 
of  ways  including 

1.  analytical  derivations  in  terms  of  obtainable  or  estimable 
inputs 

2.  simulation  (war  gaming) 

3.  extrapolation  of  empirical  data  from  operational  or 
development  tests 

k.      a  combination  of  the  above 

The  model  user  may  choose  any  of  these  methods  so  long  as  the 
resulting  inputs  are  compatible  with  the  allocation  model  „  However^, 
any  of  the  methods  used  will  require  an  assumption  about  the  tactical 
situation  with  respect  to  which  the  optimization  is  to  be  made . 
Consider  the  actual  tactical  situation  in  which  the  weapon  will  be 
used.  Assume  that  methods  exist  to  determine  the  optimum  values 
of  all  weapon  parameters,  given  that  the  weapon  will  be  used  in  a 
completely  specified  tactical  situation*  The  optimum  parameter 
values  determined  will  be  different  for  each  different  tactical 
situation  considered.  The  space  of  different  tactical  situations 
is  infinite.  However  a  single  set  of  parameter  values  must  be  spe° 
cified  for  the  weapon,  ultimately  to  be  produced  in  quantity „  By 
what  criteria  is  the  specific  tactical  situation  to  be  selected  for 
use  in  the  optimization  of  the  weapon  parameters?  The  probability 
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that  the  weapon  will  actually  be  used  in  the  tactical  situation  selected 
is  zero.   The  following  two  criteria  will  be  used  in  this  paper  in  the 
analytical  derivation  of  the  three  required  probabilities,, 
Criterion  1:   Select  an  unfavorable  tactical  situation  which  is  not 
extreme.   This  may  be  thought  of  as  a  modified  maximin  criterion., 
A  maximin  criterion  implies  optimization  of  the  weapon  parameters  to 
maximize  effectiveness  for  that  tactical  situation  which  offers 
minimum  probability  of  weapon  success.   The  extreme  situation  is 
modified  in  this  paper  because  the  low  probability  of  kill  resulting 
from  assumption  of  an  unfavorable  situation  in  the  extreme  would 
undoubtedly  distort  the  relationships  between  the  study  variables  which 
would  apply  in  a  greater  percentage  of  cases.   In  general  a  weapon 
which  is  optimized  for  an  unfavorable  situation  will  be  able  to  cope 
with  a  more  favorable  situation  for  which  its  parameter  values  were 
not  optimized. 

Criterion  2:   Select  a  tactical  situation  which  is  sufficiently 
simplified  as  to  be  within  the  scope  of  this  paper.   Progress  in 
work  of  this  kind  is  more  often  made  by  starting  with  something 
relatively  simple  and  adding  sophistication  after  experience  is 
gained  with  the  model. 


Relation  of  this  work  to  previous  work 

The  allocation  problem  has  been  studied  before  by  R„  C.  Brumfield 
in  references  1  and  2  in  which  algebraic  expressions  were  derived 
relating  total  weight  and  volume  of  a  weapon  of  arbitrary  speed  and 
range  to  performance  characteristics  and  component  characteristics „ 
This  enabled  determination  of  optimum  running  speed  and  range  for 
stern  chase  vehicles  under  conditions  of  minimum  propulsion  system 
weight;  minimum  total  weight  of  propulsion  system  and  associated 
shell;  and  minimum  propulsion  system  volume.   The  approach  in  this 
paper  is  essentially  different  from  that  of  references  1  and  2 
although  the  problem  is  much  the  same.   The  previous  work  optimized 
total  weight  and  volume  and  attack  speed  for  the  event  of  a  single 
attack  run,  given  explosive  weight,  a  length  of  run  in  the  acquisition 
phase,  and  acquisition  range.   This  paper  optimizes  an  allocation  of 
payload  weight  and  volume  to  fuel,  explosive,  and  attack  speed  for 
the  events  of  acquisition,  hit  (attack  and  reattacks)  and  kill 
given  the  maximum  allowable  payload. 

A  great  deal  of  work  has  been  done  in  the  stwdy  of  various 
phases  of  the  torpedo  problem.  References  3  and  h   form  the  background 
on  which  the  development  in  Chapter  III  is  based.   Reference  k 
develops  a  continuous  volume  scan  rate  for  a  single  circular  search 
turn.   In  this  paper,  in  the  first  three  sections  of  Chapter  III, 
a  volume  scan  for  a  single  ping  on  the  first  search  turn  is  developed. 
This  is  a  modification  of  the  continuous  scan  rate  to  the  discrete 
case  where  the  discrete  interval  is  a  single  ping  cycle.   However,  the 
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development  in  this  paper  is  continued  to  obtain  a  probability  of 
acquisition  not  only  for  a  single  circular  search  turn  but  for  con° 
tinuous  turns  until  fuel  exhaustion.,  This  further  development  is 
made  possible  because  a  tactical  situation  is  postulated  in  accordance 
with  criteria  one  and  two.  The  postulation  is  Chapter  III  limits 
the  usefulness  of  the  specifics  of  the  development  to  the  search 
pattern  considered  but  is  required  in  order  to  obtain  the  probability 
of  acquisition  input  to  the  allocation  model.  The  logic  of  the 
development  may  be  modified  as  appropriate  to  provide  for  other 
postulated  tactical  situations  in  accordance  with  requirements  to 
study  a  particular  search  pattern.  The  study  of  search  turns  subse= 
quent  to  the  first  and  until  fuel  exhaustion  is  necessary  to  provide 
probability  of  acquisition  as  a  function  of  fuel  weight  throughout 
the  range  of  the  argument . 

The  development  in  Chapter  IV  of  probability  of  a  hit  in  any  given 
pass  at  the  target  has  no  known  precedent  in  the  literature  and  is 
presented  as  an  alternative  way  of  looking  at  the  hit  probability 
problem  which  may  be  of  interest. 

It  is  emphasized  that  Chapters  III,  IV,  and  V  are  samples  of 
possible  methods  of  attacking  the  problem  of  finding  inputs  to  the 
allocation  model  of  Chapter  I  which  are  specifically  designed  to 
permit  examination  of  the  tradeoffs  between  alternative  allocations 
of  payload  weight  and  volume  among  the  three  study  variables. 
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Structure  of  the  paper 

Chapter  I  contains  the  development  of  a  mathematical  model  of 
the  torpedo  problem  in  a  form  such  that  the  dynamic  programming 
techniques  of  Chapter  II  of  reference  5  may  be  applied  to  solve  the 
allocation  decision  problem.  Chapters  III,  IV,  and  V  contain  the 
development  of  the  three  required  inputs  to  the  allocation  model  of 
Chapter  lis  probability  of  acquisition,  probability  of  hit  given 
acquisition  and  probability  of  kill  given  hit .  Appendix  A  is  a 
description  of  the  computer  flow  chart  notation  to  be  used  in  the 
flow  charts  of  the  other  appendices .  Appendix  B  is  the  computer 
flow  chart  of  the  dynamic  programming  allocation  model  of  Chapter  I. 
Appendices  C  1  through  C  8  contain  the  mathematical  details  of  the 
developments  in  sections  I  through  VIII  of  Chapter  III  plus  the 
corresponding  computer  flow  charts.  This  mathematical  supplement 
to  Chapter  in  is  included  in  appendices  to  enable  emphasis  of  the 
logic  flow  in  Chapter  III  with  a  minimum  of  mathematical  detail. 
Appendix  D  contains  the  computer  flow  chart  for  Chapter  IV. 
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Summary 

A  mathematical  model  is  provided  which  enables  optimization  of 
the  decision  problem  of  the  allocation  of  available  weight  and  space 
in  an  air°dropped  anti-submarine  torpedo  among  three  study  variables  % 
amount  of  fuel,  amount  of  explosive  and  attack  speed  capability  above 
some  preset  minimum.  The  model  is  restricted  to  the  study  of  alterna=> 
tive  allocations  which  do  not  differ  with  respect  to  their  effect  on 
weapon  reliability  and  to  their  effect  on  enemy  countermeasures 
capability  (other  than  maneuvering  capability.) 

Mathematical  models  are  developed  which  relate  (l)  probability 
of  acquisition  to  endurance  for  the  case  of  circling  search  at 
constant  depth,  (2)  probability  of  hit  given  acquisition  to  endurance 
and  attack  speed,  (3)  probability  of  hit  on  a  close~in  attack  run  to 
attack  speed,  and  (4)  probability  of  kill  damage  given  a  hit  to  amount 
of  explosive  used.  These  four  models  are  developed  as  specific 
examples  of  inputs  needed  to  be  compatible  with  the  requirements  of  the 
allocation  model.  Computer  flow  charts  are  provided  for  each  model 
except  the  last  to  facilitate  programming  in  any  desired  computer 
language  and  also  as  an  alternative  method  of  study  of  the  development 
of  the  model  logic. 
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CHAPTER  II 
THE  ALLOCATION  MODEL 
Consider  the  design  of  an  air-dropped  anti-submarine  torpedo  such 
that  the  amounts  of  fuel  and  explosive  to  be  included  and  the  capabi- 
lity  of  attack  speed  in  excess  of  some  predetermined  minimum  value 
may  be  considered  together  as  a  payload  for  which  there  is  available 
a  predetermined  maximum  of  weight  and  volume  Consider  the  employ- 
ment of  the  torpedo  from  the  time  it  splashes  into  the  water  until 
the  time  it  explodes  in  the  proximity  of  a  submarine  target  as  an 
ordered  sequence  of  the  following  three  events  each  of  which  must 
occur  in  succession  for  accomplishment  of  the  torpedo  missions 

1.  Acquisition  after  a  search  run  following  some  preset  pattern 

2.  Hit  after  an  attack  (and  if  necessary  a  series  of  reattacks) 

3.  Kill  damage  resulting  from  explosion  of  the  weapon  in  the 
proximity  ot  the  target 

Divide  the  study  variable,  amount  of  fuel,,  into  two  separate  and 
independent  study  variables: 

1.  Amount  ot  fuel  to  be  used  in  tne  event  of  acquisition  (event  1) 

2.  Amount  of  fuel  to  be  used  in  the  event  of  hit  (event  2) 
Restrict  the  study  variable,  attack  speed,  to  a  range  of  values 
such  that 

1.  The  lower  limit  of  the  range  is  some  preset  minimum  value  which 
is  greater  than  the  given  value  of  weapon  search  speed. 

2.  As  the  attack  speed  is  allowed  to  vary  over  its  range,  the  fuel 
expenditure  rate  at  attack  speed  and  at  any  depth  under  considera- 
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tion  may  also  vary.  However  as  the  propulsion  system  is  modified 
to  allow  for  changes  in  the  value  of  its  maximum  (attack)  speed 
capability,  the  fuel  expenditure  rate  at  search  speed  and  depth 
remains  approximately  constant . 
Make  the  following  assumptions; 

Assumption  1:  Weapon  reliability  and  target  capability  to  employ 
countermeasures  other  than  maneuvering  are  both  independent  of  changes 
(within  the  ranges  to  be  considered)  in  the  values  assigned  to  the 
four  study  variables: 

1.  Amount  of  fuel  for  event  1 

2.  Amount  of  fuel  for  event  2 

3.  Amount  of  explosive 
k.     Weapon  attack  speed 

Assumption  2:  The  weapon  is  100$  reliable. 

Assumption  3:   The  enemy  target  submarine  has  no  countermeasures 
capability  except  that  of  maneuvering. 

Note;  Assumptions  3  and  h   are  made  in  order  to  justify  the  relation- 
ship between  the  study  variables  and  the  measure  of  effectiveness  to 
be  introduced  later.  The  fact  that  they  are  intuitively  incorrect 
does  not  detract  from  the  validity  of  the  optimization  in  the  model 
because  of  assumption  1. 

Assumption  ki     The  three  events  of  acquisition^,  hit  and  kill  are 
independent  events  in  a  probabilistic  sense.  That  is^  if 
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P  [_xj  a  Probability  that  event  x  occurs 

P  =  P  I  target  acquisition  (event  1)1 

P„  =  P  I  target  is  hit  (event  2)  given  that  event  1  has  occurred 

P£  =  P  hcill  damage  is  inflicted  on  the  target  by  explosion  of 

weapon  warhead  (event  3)  given  that  event  2  has  occurred 

If  Pk  =  P  [target  kill] 

then  Pk  a  PA.PH-PE 

Defines 

W  =  available  pay load  weight  (lb) 

V  =  available  payload  volume  (ft  ) 

w.  =  fuel  weight  for  acquisition,  event  1  (lb) 

w1  =  fuel  weight  for  hit,  event  2  (lb) 

w  =  weight  for  modifications  to  the  propulsion  system  to  provide 

a  maximum  (attack)  speed  capability  in  excess  of  some  preset 

minimum  value  (lb) 
S  =  preset  minimum  value  of  attack  speed  (yd/sec) 
w_  =  explosive  weight  (lb) 

v.  =  fuel  volume  for  acquisition,  event  1  (ft  ) 
v'  =  fuel  volume  for  hit,  event  2  (ft  ) 
v' ' =  volume  for  modifications  to  the  propulsion  system  to  provide 

a  maximum  (attack)  speed  capability  in  excess  of  S   (ft  ) 
v_  =  explosive  volume  (ft  ) 
w  =  payload  weight  for  hit,  event  2  (lb) 


i    1 1 

w2  =  w2  +  w2 
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v2  =  payload  volume  for  hit,  event  2  (  f  t  , 


i     i  i 

v2  *  v2  +  v2 

A  w  =  increment  of  weight  of  the  any  of  the  weight  variables,  all  of 

which  are  varied  in  increments  of  equal  size  (lb) 

gi(wi)  =  the  relationship  between  the  i  th  weight  variable  and 

i   «  » 
the  i  th  volume  variable,  i  -   1,  2,  2  s    3 
t   i  i 

8i^Wi^  =  V  l  =  l>    2'    2  '  3|  81^  ■  82^ 

f,(w.)  =  natural  logarithm  (in)  of  P  expressed  as  a  function  of  w« 

f  (w' ,  w' ' )  =  In  P  expressed  as  a  function  of  w"  and  w° 9 

f_.(w0)  =  In  P_  expressed  as  a  function  of  w_, 
3  3        F  3 

3 

F  =  In  P  .    2Z  fj 

k     i  =  1  X 
i   ^  i 

W  f   W  =  lower  and  upper  limits  respectively  on  the  allowable  range 

of  fuel  weight  (lb) 
it   _i  i 
W  ,  W  =  lower  and  upper  limits  respectively  on  the  allowable  range 

of  weight  for  modifications  to  the  propulsion  system  to 

provide  a  maximum  speed  capability  in  excess  of  S  (lb) 
<  i  •   __»  «  « 

W     W 

—   '     =  lower  and  upper  limits  respectively  on  the  allowable 

range  ot  explosive  weight  (lb) 

The  concept  ot  the  lower  and  upper  limits  on  the  ranges  of  the 

study  variables  is  as  tollows.   It  will  be  obvious  to  the  design 

engineer  that  very  small  values  or  either  tuel  or  the  explosive 

weights  will  result  in  a  negligible  probability  of  success,  Thus 

there  is  no  need  to  examine  the  vicinity  of  the  extreme  values  of 

zero  and  W.   However  the  values  of  the  limits  may  be  chosen  as 

close  to  zero  and  W  as  is  desired.   In  the  case  of  W  ?    the  value 

Ik 


chosen  may  well  be  zero  as  this  implies  that  the  considered  minimum  of 

attack  speed  is  to  be  S  ,  the  preset  minimum.,   In  the  case  of  W   the 

o 

value  chosen  will  probably  be  well  under  W  as  the  maximum  of  the  range 
of  values  of  attack  speed  is  limited  by  considerations  of  resultant 
guidance  problems,  excessive  self- noise  levels,  and  so  forth.  The 
hit  event,  event  2,  contains  two  study  variables  which  must  be  con- 
sidered together.  For  events  1  and  3  only  one  study  variable  is  con- 
sidered per  event  (or  programming  stage).  Thus  for  each  value  of 
fuel  weight  (w. )  there  is  associated  a  single  value  of  acquisiti* 
probability  PA(w,  )  and  a  single  value  of  volume  v.  =  g^v^). 
Similarly  for  values  of  explosive  weight  w_. 
However 

w2  =  v'  *  V' • 


Lon 


S  » 


V2  =  V2  +  Y2     =  g2(w2}  +  g2   (W2 


This  requires  that  for  each  value  of  weight  (w  )  there  be  an 
associated  set  of  values  of  hit  probability  PH(W2^  and  of  volume 
v  .  one  value  for  each  combination  of  values  of  w  and  w   which 

Q  Q  0 

sum  to  the  considered  value  of  w  .   A  suboptimization  of  w^   and  w2 
to  provide  a  single  optimal  value  of  PH(*0  for  each  \j^   as  an  input 
to  the  allocation  model  would  neglect  the  associated  volume 
restraints.  Therefore  index  the  weight  w  by  k(k  *  1,  •••,  K.)  over 
its  range 

min  £w«  ,  W'  ^   ^ v2  ^min  £w«  +  W«  •  ,  w} 
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and  define 


w  ,  =  k  the  value  of  weight  to  be  distributed  between  w.  and 
2k  d 


wj/;  k  =!,-«.,  K  (lb) 


w'  ,  =  j  th  value  of  w'  associated  with  wfc|  j  -  1/°",  Jfc  (lb) 

w''   =  j  th  value  of  w'1  associated  with  w  ^    j  ■  1,  •••J|E  (lb) 

P    =  j  th  value  of  P  associated  with  w  ,  j  j  -  1/",  J, 

v'  ,  =  j  th  value  of  v'  associated  with  w  ,  j  j  =  ls    °°°,    J^°s 

k  =  1,  •-.,  K 
v''   =  j  the  value  of  v'1  associated  with  w  ,  j  j  -  1,'»«,  Jfe| 

v  ,  =  v'  ,  +  v'1.,  =  j  the  value  of  v^  associated  with 
2jk    2jk    2jk   J  2 

w2k;  j  =  1,  ••,  Jk  (ft3) 

The  available  payload  W  and  the  upper  and  lower  weight  limits  or 
the  weight  increment  A  w  are  so  chosen  that  all  the  weights  are 
divisible  by  Aw.  Then  there  will  be 


Aw     Aw 


considered  values  of  w,  denoted  w.  j  r  ■  1,  ' ' ' t   R, 
there  will  be 


min  4  W'  +  W» ' ,  W  j-    min  £  W»  9  W'  *j- 
Aw  Aw 


=  K 


considered  values  of  w  denoted  w  ,  ;  k  =  1,  °°°?   K 
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and  there  will  be 

Aw     Aw 
considered  values  of  w^  denoted  w   j  s  =  1,  O0%  S 

j  jS 

The  functions  g1(w1),  g^Cw^),  g^'  (w^'  ),  83^3)*  f^Wj),  and 
f  (w^)  are  assumed  to  be  non-decreasing .  That  is,  an  increase  in 
the  argument  will  not  decrease  the  value  of  the  function „   In  the 
case  of  f  (w  ,  w   )  it  is  assumed  that  any  increase  in  either  of 
the  arguments  while  maintaining  the  other  fixed  will  not  decrease 
the  value  of  the  function. 

The  problem  is  to  maximize  F  =  In  P, 
subject  to 

vl  +  w2  *  W3  -W 
vl  +  v2  *  v3  ~V 

Consider  the  new  problem  to  maximize 

F»  =  F  -  X  [g1(w1)  ♦  gj(wj)  ♦  g^  (wj«)  ♦  g3(w3)J 


subject  to 


w.  +  w  +  w_  <W 


gi(wi)  ^0;  i  =  1,  29  2  ,  3 
where  X  is  a  non-negative  real  number  (Lagrange-multiplier ) 

Define  the  recurrence  relations' 

hx(w)  =  fx(w)   -  X  gl(w)         if  W  =  w^; 

=  -  oo  otherwise 
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r  = 


h2(W)  = 


max  f"   max        [f   (w*         v'«    )  -  X  w  .J] 

w     ;  k  =  1,   •••,  K  1  1  <j  <J,    L  z     ZJJC      ZJJ<:  23K2 


where 


f2Hlk'  "Ilk*  "  *  PH1k  ^  *  FH(w2ik>  W2ik} 


2v"2jk9  "2jk 


Hjk 


V2jk=S2(w^k)*gJ'(wJjk) 


h3(W)  = 


max 
w3s»  s  =  1>  *"'  S 


f3^w3s^  "  Xg3^w3s^  +  V 


=V] 


The  dynamic  programming  methods  of  Chapter  II  of  reference  5  are 
next  applied  to  the  recurrence  relations  h.,  h  ,  and  h  for  some  fixed 

X  to  obtain  a  solution  for  the  optimum  value  of  h  (w)  ■  F  and  for  the 
corresponding  optimum  values  of  w. ,  w  ,  w   ,  and  w_  . 

Next  the  corresponding  optimum  values  of  v  .,  v  -  and  v  are  com= 
puted  using  the  computer  determined  optimum  values  of  w. ,  w  ,  w  s    and  w  , 
The  sum  (v^  +  v  +  v  )  is  then  compared  with  the  available  payload 
volume  V.   Probably  v..  +  v  +  v  f   V. 
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Then  if 

vl  +  v2  +  v3>  V 
the  value  of  X  must  be  increased  and  the  dynamic  program  repeated 
until  the  minimum  value  of  X  is  found  which  results  in 

Vl  +  V2  +  V3  ~V 
If  v.  +  v  +  v  <V 

the  value  of  X  is  decreased  until  the  minimum  value  of  X  is  found 

which  results  in 

Vl  +  V2  +  v3  ~V 

The  optimum  value  of  h  (w),  w, ,  w  ,  w'1  and  w_  determined  for 

such  a  X  specify  the  solution  values  of  the  three  study  variables . 

The  Lagrange-multiplier,   X  >   may  be  thought  of  as  a  cost 

coefficient.   If:   X=  ° 

Vj_  +  v2  +  v^  <V 

then  the  volume  restraint  may  be  thought  of  as  having  cost  nothing . 
That  is,  the  measure  of  effectiveness  has  been  optimized  solely  with 
respect  to  the  weight  restraint  and  the  result  has  met  the  volume 
requirements.  On  the  other  hand,  if  X-   o 

Vj  +  v  +  v_>  V 
then  some  positive  value  must  be  assigned  to  X  which  will  result 
in  meeting  the  problem  volume  restraint  v^  +  v  +  v  <V„ 

Ultimately  the  subsequent  process  of  optimization  with  some 
positive  value  of  X  will  result  in  smaller  optimum  values  of  one 
or  more  of  the  weight  variables  with  a  correspondingly  reduced  value 
of  the  measure  of  effectiveness. 
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This  reduction  in  the  solution  value  of  the  measure  of  effective^ 
ness  may  be  considered  as  the  cost  of  having  a  volume  restraint „  Of 
course  if  weights  and  volumes  are  interchanged  in  the  model  the 
resulting  X  will  be  a  measure  of  the  cost  of  the  weight  restraint „ 
If  the  solution  value  of  probability  of  kill  is  desired  then 
Pk  =  exp  [h3(w)  +  X(vx  +  v2  +  v3)j 
where  h  (w),  v.,  v  and  v  are  solution  values. 

The  computer  flow  chart  to  accomplish  this  dynamic  programming 
procedure  is  contained  in  Appendix  B. 

It  may  be  of  interest  to  allow  the  values  of  the  available 
pay load  weight  and  volume  (W  and  V)  to  vary  in  the  vicinity  of  their 
preset  values.   Then  the  sensitivity  of  the  effectiveness  measure^, 
P.  ,  to  such  small  variations  may  be  examined  in  order  to  determine 
the  possible  value  of  modifications  to  the  weapon  delivery  systems  to 
allow  for  a  different  (more  effective)  load  capability „ 

Input  Requirements  for  the  allocation  model 

W  =  Available  payload  weight  (lb) 

V  =  Available  payload  volume  (ft  ) 

W  ,  W  =  lower  and  upper  limits  respectively  on  the  allowable  range 

of  fuel  weight  (lb) 
ii   _n 
W  ,  W   =  lower  and  upper  limits  respectively  on  the  allowable 

range  of  weight  for  modifications  to  the  propulsion 

system  to  provide  a  maximum  (attack)  speed  capability  in 

excess  of  a  minimum  preset  value  (lb) 
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Ill   _lll 
W   ,  W    =  lower  and  upper  limits  respectively  on  the  allowable 

range  of  explosive  weight  (lb) 
1 1 
g,(w.);  i  =  1,  2  ,  3  =  the  relationships  between  the  i  th  weight 

variable  and  the  i  th  volume  variablej 

1 1 
i  =  1,  2  ,  3 

P  (w. )  =  probability  of  target  acquisition  expressed  as  a  function 

of  the  fuel  weight  for  acquisition 

P  (w  ,  w   )  =  probability  of  hitting  the  target  given  that  target 

acquisition  has  occurred  expressed  as  a  function  of  the 

fuel  weight  for  hit  and  of  the  weight  for  modifications 

to  the  propulsion  system  to  provide  a  maximum  (attack) 

speed  capability  in  excess  of  S 

P  (w  )  =  probability  of  target  kill  given  that  a  hit  has  occurred 

expressed  as  a  function  of  the  explosive  weight 

S  =  minimum  preset  value  of  attack  speed  (yd/sec) 

Aw  =  size  of  the  increments  of  weight  of  the  weight  variables^ 

all  of  which  are  varied  in  increments  of  equal  size  (lb) 

The  payload  weight  and  volume  are  determined  by  the  capabilities 

of  the  weapon  delivery  systems.  The  limits  of  the  weight  variables 

are  determined  by  a  preliminary  study  which  reflects  limiting  design 

considerations  and  reduces  the  possible  ranges  of  the  weight 

variables  to  reduce  computer  time«  The  weight=volume  relationships 

are  determined  by  engineering  considerations  reflecting  the  state 

of  the  art.  The  minimum  attack  speed  is  determined  by  intelligence 

estimates  or  CNO  requirements.  The  size  of  the  weight  increments 
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is  determined  by  considerations  of  required  accuracy  of  results  and 
available  computer  time.  The  probabilities  of  acquisition,  hit  and 
kill  may  be  obtained  in  a  number  of  ways  including 

1.  analytical  derivations  in  terms  of  obtainable  or  estimable  inputs 

2.  simulation  (war  gaming) 

3.  extrapolation  of  empirical  data  from  operational  or  development 
tests 

k.     a  combination  of  the  above 

The  remaining  chapters  provide  samples  of  possible  analytical  methods 
of  providing  the  three  required  probability  inputs  which  are  specifi= 
cally  designed  to  be  compatible  with  the  allocation  model  requirements 
It  is  possible  to  extend  the  model  to  provide  for  cost  restraints 
The  procedure  would  require  a  second  Lagrange=multiplier  in  the 
function  to  be  maximized  (F  ).  Additional  inputs  required  would  be 
the  relationsuips  between  the  costs  and  the  corresponding  weight 
variables.   It  is  considered,  however,  that  cost  restraints  are  not 
significant  to  the  problem  considered  in  this  paper „  Variations  in 
the  amount  of  fuel  and  amount  of  explosive  have  little  effect  on  the 
costs  involved.  Even  changes  in  attack  speed  within  the  limits 
appropriate  for  consideration  will  probably  have  negligible  cost 
effect  when  compared  to  the  total  cost  of  developing,  producing  and 
operating  a  weapon  in  the  course  of  its  life  cycle. 
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CHAPTER  III 
PROBABILITY  OF  ACQUISITION 
Consider  a  tactical  situation  in  which  a  weapon  is  delivered  by 
some  air  delivery  weapons  system  to  a  splash  point  on  the  surface  of 
the  water  in  the  vicinity  of  a  submarine  target  and  make  the  following 
assumptions , 

Assumption  1;   The  splash  point  defines  the  center  of  the  top  of  a 
right  circular  cylinder  called  the  target  volume  cylinder .  The 
probability  is  quite  small  and  is  assumed  zero  that  the  target  is  out- 
side of  the  target  volume  cylinder.  The  target  is  placed  at  random 
within  the  target  volume  cylinder. 

Assumption  2:  The  weapon  dives  to  its  assigned  search  depth  on  a 
straight  line  track  and  commences  a  right  circular  search  turn  at 
constant  depth,  speed,  and  turn  rate  immediately  upon  reaching 
search  depth.  The  weapon  searches  in  the  active  acoustic  mode 
commencing  its  first  ping  upon  commencement  of  the  search  turn. 
The  weapon  continues  its  turn  until  target  acquisition  or  fuel  ex- 
haustion. 

Assumption  3:   The  weapon  does  not  acquire  a  false  target. 
Assumption  k:      The  search  turn  of  the  weapon  is  approximated  by 
rotating  the  weapon  about  a  vertical  axis  through  the  point  at  which 
the  weapon  commences  search. 

Assumption  5:  While  the  weapon  is  diving  to  search  depth,  the  radius 
of  the  target  volume  cylinder  increases  at  a  rate  equal  to  the  target 
cruise  speed. 
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Assumption  6;  The  target  is  alerted  on  commencement  of  the  weapon  search 

turn  by  the  first  'ping'  of  the  weapon.  Thereafter  the  target  opens  at 

constant  depth  and  speed  radially  from  the  projection  of  weapon  position 

on  the  horizontal  plane  through  the  target. 

Assumption  7;   During  search  the  radius  of  the  target  volume  cylinder 

increases  at  a  constant  rate  equal  to  the  target  speed. 

Assumption  8;   The  'definite  range'  law  as  developed  in  Chapter  II  of 

reference  6  applies  to  acquisition.  Thus  on  a  given  radial  during  a 

ping  cycle  under  specified  conditions  there  is  a  range  beyond  which 

acquisition  probability  is  zero  and  within  which  acquisition  probability 

is  one.  Regions  closer  than  the  'definite  range"  limits  are  considered 

scanned , 

Assumption  3%      During  search  the  weapov  makes  a  series  of  independent 
glimpses  for  the  target,  one  for  each  ping  cycle. 

Assumption  10;  Target  cylinder  volume  is  considered  constant  during 
each  ping  cycle  at  the  value  valid  at  the  initiation  of  the  ping. 
Assumption  11:  The  water  medium  is  considered  to  be  isotropic. 
That  is,  all  environmental  parameters  that  affect  acoustic  performance 
have  a  constant  value  throughout  the  considered  volume. 
Assumption  13;  Consider  planes  containing  the  transducer  and  angles 
measured  at  the  transducer.  Consider  the  vertical  plane  containing  the 
weapon  longitudinal  axis  as  the  vertical  reference  plane.  As  the 
weapon  sweeps  through  an  angle  0  in  the  horizontal  plane  it  sweeps 
through  a  lesser  angle  0t  in  a  plane  perpendicular  to  the  vertical 
reference  plane  but  inclined  to  the  horizontal.  Also  if  the  inclination 
of  the  plane  inclined  to  the  horizontal  as  measured  in  the  vertical 
reference  plane  is  «   then  the  inclination  of  the  plane  inclined  to  the 
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horizontal  as  measured  in  a  vertical  plane  other  than  the  reference 
plane  is  a  lesser  angle  «' .   In  this  paper  it  is  assumed  that 

0«  =  0 

ec«  =  ex 
In  practical  situations  likely  to  be  considered,  this  assumption 
results  in  much  smaller  errors  than  those  created  by  treating  the 
range  problem  in  discrete  steps  as  is  done  in  this  chapter. 

Define : 

P  =  P  (i)  =  P  T acquisition  in  I  ping  cycles  J 

where  P  I"  x]  =  probability  of  event  X  happening 

g.  =  P  [acquisition  during  the  i  th  ping  cycle J  $  i  =  1, 

w   =  weapon  fuel  expenditure  rate  at  search  depth  and  speed 

PRR  =  weapon  active  sonar  pulse  repetition  rate  for  search 
s 

(pings/sec) 

w    =  fuel  used  by  weapon  to  reach  search  depth .  (lb) 
OSD 

w  =  fuel  available  for  acquisition,  (lb) 
Then,  as  developed  in  Chapter  2  of  reference  6, 
I 

pA  =  i  -  ~TT  d-gi) 

A       1=1 


sec 


and  since 


<wl  ~  W  PRRs 


'It 


=  I 


and  '[jc]  =  greatest  integer  ^X   and  where  I 
is  the  number  of  complete  ping  cycles  available  for  the  event  of  acquis!* 
tion,  it  may  be  said  that  P.  is  a  function  of  W]L  and  a  table  may  be 
computed  of  PA  vs  Wj  where  Wj  is  allowed  to  vary  between  W  and  W  . 
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Of  course  if  I  is  less  than  1  then  P  is  defined  to  be  zero.  The 

probability  of  acquisition  in  a  specified  glimpse  (gi)  may  be  con= 

sidered  as  a  coverage  factor.,  The  volume  containing  the  target  with 

positive  probability  is  a  region  to  be  covered,  A  portion  of  the 

region  to  be  covered  is  scanned  (covered)  by  the  weapon  in  accordance 

with  the  definite  range  law.   Because  the  target  has  been  randomly 

placed  in  the  target  volume  cylinder  (region  of  target  uncertainty 

to  be  covered)  then 

volume  scanned  on  i  th  ping  cycle 
8i  "  target  volume  during  i  the  ping  cycle 

Thus  the  problem  is  reduced  to  finding  what  portion  of  the  volume 
of  target  uncertainty  is  scanned  by  the  weapon  on  any  ping.  This  will 
be  different  for  every  ping. 

The  following  development  for  P  is  divided  into  eight  consecutive 
steps.  The  development  is  amplified  as  appropriate  and  complete  flow 
charts  for  a  computer  program  are  provided  in  appendices  C  1  through  C  I 
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1.  Available  Echo  Excess  (B,k) 

Consider  a  vertical  plane  containing  the  longitudinal  axis  of  the 
weapon  at  the  time  of  ping  transmission .  Consider,  as  in  Figure  1,  a 
second  plane  perpendicular  to  the  first  and  containing  the  weapon 
transducer  such  that  it  makes  an  angle  j  with  the  weapon  longitudinal 
axis,  A  set  of  such  planes  is  considered  such  that  the  angle  j  is 
associated  with  the  j  th  plane .   In  the  vertical  plane  angles  are 
measured  positive  counterclockwise  from  the  reference  axis» 

Design  or  operational  considerations  will  make  it  necessary  to 
prohibit  reception  of  signals  from  off  the  transducer  axis  in  the 
vertical  plane  in  excess  of  some  limting  angle  J,  not  greater  than 
90  degrees.  Allowing  one  degree  between  each  of  the  considered  planes  j, 
the  index  j  will  assume  the  range  of  values  -J,  -J  +  1*  •••#-!*  °»    l» 
...,  J.  Thus  2J+1  planes  are  to  be  considered,  one  containing  the  trans- 
ducer  axis  and  J  planes  on  either  side. 

Consider  increments  of  target  range  of  100  yards  in  the  plane  j. 
For  an  assumed  target  positioned  in  the  j  th  plane  at  a  range  of  100k 
yards  on  the  transducer  azimuth  axis  projection  in  the  plane  j  (point  P 
in  Figure  l),  there  exists  an  echo  excess  (Bjk)  by  which  the  signal  level 

(SL   )  of  the  echo  reflected  from  the  target  exceeds  the  background 

v   jk' 

masking  level  at  the  weapon  transducer  (M(BG)k).   That  is i 

B.,  =  SL..-  M(BG).  J  j=  -J,  "•»  J°>    k=1>  '">    KJ 

jk       JK.         K 

Assuming  a  definite  range  law  of  acquisition,  a  negative  echo  excess 
yields  zero  probability  of  acquisition  and  a  positive  echo  excess  yields 
a  probability  of  acquisition  of  one.  The  maximum  range  of  acquisition  on 
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Figure  lj  The  plane  j.  The  transducer  is  at  the  origin  pointing  in 

x  direction.  Line  ox  is  the  transducer  axis  extension .   Lines 
ox,  oP,  and  oy  are  in  the  vertical  plane.   Lines  oP  and  oz  are 
in  the  plane  j.  The  polar  coordinates  of  the  point  P  in  the 
vertical  plane  using  the  x  axis  as  a  reference  are  (100k,  j)o 
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the  azimuth  axis  in  the  j  th  plane  is  thus  the  range  of  100K.  yards  at 

which  B   equals  zero.  At  ranges  less  than  100K8  yards  there  is  an 

jk  .         J 

echo  excess  of  B.,  above  that  required  for  acquisition  on  the  azimuth 

axis.  This  excess  is  available  for  use  off  the  aximuth  axis . 

The  signal  level  varies  with  assumed  target  range  and  with  the 
elevation  angle  of  the  target  above  or  below  the  tranducer  axis  in 
accordance  with 

SL.k  =  WS  +  TS  -  2  Hfe  ♦  10  log  [^(J)  b2(j)J 

where: 

WS  =  on-axis  source  level  (db  at  1  yard  from  transducer) 
TS  =  target  strength  of  target  (db) 

H  =  one  way  transmission  loss  at  a  range  of  100k  yards  (db) 

k 
b  (j)  =  normalized  vertical  pattern  function  for  transmitting 

pattern  (i.e.,  the  ratio  of  intensity  at  the  angle  j 
off  the  vertical  axis  but  on  the  azimuth  axis  to  intensity 
on-axis  in  both  the  vertical  and  aximuth  planes, 
b  (j)  =  normalized  vertical  pattern  function  for  receiving  pattern. 
There  are  several  different  kinds  of  masking  noises  and  it  is  their 
total  effect  which  is  considered  as  the  interference  background  masking 
level.   Included  are  self-noise,  target  noise,  ambient  noise,  and 
volume  and  boundary  reverberation. 

There  is  also  an  enabling  level  designed  into  the  weapon  electrical 
circuitry  to  protect  the  acoustic  system  from  excessive  interference 
and  to  avoid  interpretation  of  the  interference  as  a  signal.  As  used 
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here,  the  weapon  enabling  level  (WEI^)  includes  the  weapon  recognition 
differential,   It  is  a  function  of  range  (or  time)  and  decreases  with 
increasing  range.  The  background  masking  level  (M(BG)k)  is  the  larger 
of  the  two:   interference  background  masking  level  and  the  weapon 

enabling  level. 

The  output  of  this  step  is  a  table  of  the  available  echo  excess 

(s.k;  j=  -J,   •••*  J;   k=  l>   °'°>   Kj)' 
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2.  Scanned  Volume  Radials;  General  Case 

In  each  of  the  planes  j  there  is  now  defined  for  each  value  of 
range  (100k)  an  available  echo  excess  Bjk  for  acquisition  on  the 
azimuth  axis.  This  available  excess  may  be  used  to  find  limits  for 
off-axis  acquisition.  Define  an  angle*  fc  to  be  the  angle  through 
which  the  weapon  turns  while  a  ping  travels  to  a  range  of  100k  and 


returns. 
Then 


200ko}g 
k  ~~  c 

where 


»k- 


CD 
S 


=  weapon  search  turn  rate  (deg/sec) 
c  =  velocity  of  sound  in  the  medium  (ft/sec) 
For  echo  ranging  in  the  plane  j  to  a  range  of  100k,  an  energy  pulse 
transmitted  at  angle  p  will  return  at  the  angle  p  -0fc. 
The  echo  excess  A.,  (p)  required  to  acquire  the  target  at  range  100k 
in  plane  j  at  angle   p  off  the  aximuth  axis  is 

A  (  p)  =  -  10  log  b3  ( P)  -  10  logb^(  P  -*  k) 
where  b  (  p)  -  normalized  aximuth  pattern  function  for  transmitting 
pattern 
b,  (  P)  =  normalized  aximuth  pattern  function  for  receiving 

pattern 
If  the  required  echo  excess  A jk  (  p)  exceeds  the  available  echo  excess 

B   then  acquisition  will  not  take  place  at  range  100k  at  an  azimuth 

Jk 
removed  p  degrees  from  the  reference  axis  of  plane  j.  However  if 

A   (  B)  does  not  exceed  B,.,  acquisition  takes  pl< 
jkN  r  J* 
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Lace, 


In  each  plane  j  and  at  each  target  range  100k  there  are  two  limiting 

angles  between  which  acquisition  will  occur  with  probability  one  and 

outside  of  which  there  will  be  no  acquisition .   Looking  from  the  trans- 

ducer  toward  the  potential  target  define  the  set  of  limiting  angles 

B    with  m  =  1  for  the  limiting  angle  to  port  and  m  ■  2  for  the 
Kjmk 

limiting  angle  to  starboard.  As  shown  in  Figure  2,  the  reference 
axis  for  polar  coordinates  in  the  j  th  plane  is  the  projection  of 
the  transducer  axis  at  the  time  of  ping  transmission  in  that  plane . 
Positive  angles  are  measured  clockwise  from  the  reference  axis  in  all 
planes  j  and  in  any  horizontal  plane.  The  pattern  may  not  be  centered 
on  the  axis  and  thusp\mk  may  be  positive  or  negative.  Design  or 
operational  considerations  will  make  it  necessary  to  prohibit  reception 
of  signals  from  off  the  transducer  axis  in  the  azimuth  plane  in  excess 
of  some  limiting  angle  «  in  the  direction  opposite  to  the  search  turn 
and  <r  '  in  the  direction  of  the  search  turn.  This  restricts  the 
angles  p\  .  as  follows 

B    -  0   _  « 
Pjlk   \\ 

I  ft     -<t     ^oc  ' 
I  Pj2k   *  k  | 

Recall  that  in  plane  j  the  maximum  range  for  a  non-negative  B^k 
was  computed  to  be  100K1.  yards.  For  a  straight  running  weapon  this 
would  be  the  maximum  acquisition  range  in  plane  j.  Actually  an 
acquisition  range  of  100K1.  yards  will  not  be  realized  in  plane  j  in 
this  problem  because  the  weapon  is  turning  at  some  rate.  The  energy 
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from  a  given  source  direction  will  never  return  as  an  echo  in  the  same 
direction.  Thus  there  will  always  be  some  penalty  for  off-axis  trans- 
mission or  reception. 

Unless  the  turn  rate  is  extremely  high,  the  insonified  volume  of 
adjacent  weapon  pings  will  overlap.   If  a  given  ping  scans  part  of  the 
volume  insonified  by  the  previous  ping,  the  overlap  cannot  be  credited 
to  the  volume  scanned  by  the  later  ping  because  the  overlap  volume  was 
removed  from  the  volume  of  target  uncertainty  by  the  previous  pingo 
The  concept  of  increased  probability  of  detection  because  of  overlap 
is  not  valid  in  this  case  because  of  the  use  of  the  definite  range 
law.   Thus  all  overlap  volume  should  be  credited  to  the  ping  which 
first  scanned  it.   This  would  give  all  pings  subsequent  to  the  first 
ping  during  the  first  search  turn  an  equal  scanned  volume  less  than 
that  of  the  first  ping.   In  this  paper  the  situation  is  approximated 
by  dividing  the  volume  of  overlap  between  the  pings  concerned  to  give 
all  first  turn  pings  an  equal  volume  of  insonif ication.  This  approxima= 
tion  decreases  P  (l)  but  the  error  in  P  (i)  decreases  with  increasing 
I  until  the  end  of  the  first  search  turn  when  it  becomes  negligible  „ 

The  plane  of  maximum  overlap  is  plane  0.   In  this  plane^  for  small 
k,  Pnp.  for  ping  i  will  be  larger  than  b  +  p      for  ping  (i+l). 
See  Figure  3«  Here  o  is  the  angle  between  the  axes  of  pings  i  and 
(i+l)  or  the  angle  through  which  the  weapon  turns  between  pings „  As  k 
becomes  larger  the  differences  between  P   ,  for  ping  i  andb  +  p_,,  for 
ping  (i+l)  become  smaller.   Define   p1  to  be  the  angle  at  which  the 
two  are  equal.  Thus  the  right-most  angle  for  acquisition  in  any  plane 

3^ 


(1*1) 


Figure  3:   Overlap  of  insonified  volume  between  adjacent  pings 
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is  not  greater  than  p1  and  the  leftmost  angle  is  not  less  than  p!  -  & . 

Some  of  the  extreme  planes  j  may  not  be  useful  for  acquisition  at  all. 

By  including  for  consideration  only  those  planes  in  which  acquisition 

may  take  place  the  index  j  takes  on  the  new  range  from  J1  below  the 

axis  to  j'  '  above  the  axis . 

The  range  or  radial  distance  associated  with  the  angle  pjmk  is 
100k  yards.   See  Figure  k.     There  is  a  range  maximum  called  the  ping 
interval  range  defined  by  the  time  between  pings  (l/PRR8)  and  the 
velocity  of  sound  in  the  medium  (c).  Further,  for  planes  inclined 
above  the  horizontal  there  is  another  range  maximum  defined  by  the 
range  at  which  the  radial  at  angle  ^   breaks  the  surface  of  the  water 
Similarly  for  planes  inclined  below  the  horizontal  there  is  a  range 
maximum  defined  by  the  range  at  which  the  radial  at  angle  pjmk  breaks 
the  maximum  target  depth  (D^).   Volume  scanned  below  the  maximum 
target  depth  does  not  contribute  to  acquisition  probability.  Thus  the 
acquisition  range  at  any  angle  Pjmk  is  the  minimum  of  the  ranges  just 

discussed  to  be  denoted  Rjmk« 

Since  the  index  k  is  no  longer  of  value  to  indicate  range  it  is 
useful  now  to  relabel  the  radials  such  that  the  index  k  indicates  the 
separation  from  the  leftmost  radial  in  the  plane.  See  Figure  5.  The 
leftmost  radial  is  labeled  (RjL,  Pjl).  The  general  radial  is 

(R,,  PJ.  The  rightmost  radial  is  (R   ,  p„).     For  planes  close 
v  jk'  rjk  J  j     j 

to  the  axis,  probably 

Vj  -  » '   and   V. -  'n  °b 

This  may  not  be  true  of  the  extreme  planes. 
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Figure  k:       Acquisition  range  at  angle   p  fe.   Line  OR  has  length 
100k.  Line  SS1  denotes  the  surface.   Length  01  is  the 
ping  interval  range.  Lined  area  is  scanned. 
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jypjKj 


(Rjk,pjk) 


Figure  5:   Relabeled  radials  in  plane  j 


38 


The  output  of  this  step  is  tables  of 
Pjk  and  R.ki  j  =  J1,   •',  JMi  k  =  L>  °'°>   Kj 
plus  the  angles  p'  and  O  . 
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3.   Scanned  Volume:  General  Case 

The  radials  computed  in  step  2  define  areas  of  insonification  in 
each  of  a  set  of  planes  j.  The  area  in  each  plane  may  be  computed  as 
the  sum  of  a  set  of  triangles  as  shown  in  Figure  6. 

y1 

a  =  (1/2)  S=  Rjk  RJ(k+i)  sln  fy(k*i)  -  Pjk 

J         k=l 


.) 


The  irregular  shape  of  the  areas  in  the  planes  j  makes  it  difficult 
to  work  with  them  in  their  present  form.  Define  a  set  of  isosceles 
triangles  (one  for  each  plane  j)  such  that  the  triangle  areas  are 
equal  to  the  scanned  areas  A.  in  planes  j  and  the  common  apex  angle 
is  the  angle  turned  by  the  weapon  per  ping  (b).   See  Figure  ?.  The 

altitude  of  each  triangle  is 

11/2 


*r 


AJ 


tan(b/2) 


The  set  of  triangles  so  defined  approximates  the  areas  scanned  and 
makes  it  possible  to  obtain  the  total  scanned  volume  SVQl  on  ping  zero 
on  the  first  turn.   Ping  zero  is  a  ficticious  ping  denoting  the 
general  case.  The  computation  of  SV±l   (scanned  volume  on  i  th  ping  of 
first  search  turn)  is  identical  to  that  for  SVQl  to  be  developed  next. 

Define:  |i  =  weapon  angle  of  attack  (pitch)  (deg). 
The  plane  j  makes  an  angle  of  (j+ji)  degrees  with  the  horizontal. 

The  set  of  isosceles  triangles  is  used  to  define  a  set  of  pyramids 
A  vertical  cross-section  of  the  set  of  pyramids  is  shown  in  Figure  8 


>(Rjk»Pjk) 


Rj(k+ir  Pj(m)  } 


Figure  6:   Area  of  insonif ication  in  plane  j 


kl 


Figure  7:   Equivalence  of  area  A.  and  isosceles  triangle 
OHH1  of  altitude  H.  and  apex  angle  b  . 
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Figure  8:   Vertical  cross~section  of  insonified  volume  on 
ping  zero,  x  is  the  horizontal  axis,  y  is  the 
vertical  axis,  and  OQ  is  the  transducer  axis 
extension. 
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and  a  single  pyramid  is  shown  in  Figure  9°  The  volume  of  a  single 
pyramid  is 

v.  =  l/3    area  of  base  |  J  altitude! 

As  shown  in  Figure  9>  2X  is  the  average  value  of  one  of  the 
base  dimensions  and  Y  is  the  value  of  the  other  dimension. 
Thus  the  base  area  is  2XY. 


X  =  [(Hj  ♦  H(j+a))  tan(6/2)]/2 

0 

Since   sin0  =  h/fe.  =  (H/^n  sin  1  )/l 


the  altitude  is 

h  =  (H.  H(j+1)  sin  1°)A 

and 

7.  =  (l/3)(2XY)h  =  (V3)(Hj4Hy+1))  Hj  H(j+l)  sin  1°  tan(&/2) 

The  volume  scanned  is 


SV01=  S=  Vj 


The  output  of  this  step  is  tables  of  H.£  j  =  J  s    ■ • »,  J 


i  0 


and  the  volume  SV  .  . 
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k .   Special  Volume  Elements;  General  Case 

For  the  case  of  search  at  constant  depth  in  a  continuous  circle 
as  is  being  considered  in  this  paper,  the  situation  changes  drastically 
upon  completion  of  the  first  search  turn  circle „   Subsequent  search 
turns  will  effectively  scan  much  less  volume  because  the  weapon  retraces 
its  path  of  the  previous  turn.  Define  the  maximum  available  number  of 
search  turns  (u)  to  be 


U  = 


*  -  W0SD)  Ms^  Wlt] 


+  1 


where  M  is  the  number  of  ping  cycles  in  one  search  turn 
M  =   II360/&II   +  1 


H' 


greatest  integer  ^  x 


=   greatest  integer  <  x 
W   =  maximum  allowable  weapon  fuel  (lb) 

w.   =  fuel  expenditure  rate  at  search  depth  and  speed  (lb/sec) 
w    =  fuel  used  in  dive  to  search  depth  (lb) 


(^  Dws  »ot 


w 


080         sin0DS  SWD 


l/3  D   =  weapon  search  depth  (yd) 

W  o 

w   =  average  weapon  fuel  expenditure  rate  in  dive  to  search 

depth  (lb/sec) 
S^  =  weapon  speed  in  dive  to  search  depth  (yd/sec) 
0    =  angle  of  dive  to  search  depth  (deg) 
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Consider  a  side  view  or  vertical  cross-section  of  the  scanned  volume 
of  a  single  ping  as  in  Figure  10.  The  vertically  lined  area  is  that 
scanned  by  the  weapon  on  the  i  th  ping  of  the  first  search  turn. 
Recall  assumption  6  that  the  target  opens  radially  from  the  weapon 
position  and  assumption  12  about  weapon  orientations  on  different 
search  turns.  Define  a  distance  of  TES  yards  as  the  distance  traveled 
by  the  target  during  the  time  it  takes  the  weapon  to  complete  one 
search  turn.  The  weapon  scanned  area  on  the  i  th  ping  of  the  second 
search  turn  is  equivalent  to  that  of  the  areas  lined  sloping  downward 
to  the  right  on  Figure  10.  The  weapon  actually  scans  the  same 
geographic  volume  on  each  search  turn  but  since  the  target  is  assumed 
to  open  radially  the  situation  is  equivalent  to  translating  the 
area  scanned  on  the  first  turn  a  distance  TES  to  the  left  for  each 
subsequent  turn.   If  the  target  is  in  the  area  lined  sloping  downward 

o  n 

to  the  left  on  the  i  th  ping  of  the  first  search  turn^  he  will  escape 
prior  to  the  i  th  ping  of  the  second  search  turn.   Similarly  for  sub= 
sequent  search  turns.  A  target  above  the  critical  depth  CDA  or  below 
the  critical  depth  CDB  may  escape  but  a  target  at  any  intermediate  depth 
will  be  acquired  eventually  if  fuel  supply  permits  because  of  the  assump= 
tion  of  an  opening  target.  The  search  turns  after  the  first  are  desirable 
because  of  the  possibility  of  acquiring  a  target   trapped   inside  the 
insonified  volume  of  the  first  search  turn.  Treatment  of  scanned 
volumes  on  the  second  and  subsequent  turns  will  differ  depending  on 
whether  the  volume  is  above  or  below  the  depths  CDA  and  CDB.   Although 
the  areas  scanned  above  depth  CDA  on  search  turns  1  and  2  are  equal  in 
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Figure  10,  the  volumes  represented  by  these  areas  are  not  equal .  A  top 

view  of  ping  i  is  a  pie  shaped  slice  which  makes  the  volume  associated 

with  the  first  search  turn  area  larger  than  that  of  the  corresponding 

second  search  turn  area.  Thus  each  volume  element  must  be  calculated 

independently  and  the  volume  scanned  by  the  weapon  on  ping  i  of  search 

turn  u  (SV.  )  is  the  sum  of  several  independently  calculated  volumes . 

In  this  step  the  volume  elements  in  the  general  case  (i=0)  are  computed „ 

These  are  modified  as  necessary  for  each  ping  orientation  in  a  later 

step.  Although  Figures  11  and  12  and  the  development  to  follow  involve 

the  volumes  above  the  search  depth,  the  development  for  volumes  below  the 

search  depth  is  identical. 

Consider  the  vertical  cross-section  of  Figure  10.  The  cross= 

section  has  been  split  into  several  regions  or  volume  layers  identified 

on  the  left  margin  by  their  bottom  boundary  labeled  L,  such  that  the 

bottom  boundary  of  a  level  (or  layer)  at  or  above  depth  CDA  has  a 

length  of  L,  inside  the  first  turn  scanned  region.  The  planes  j  are 
d 

considered  using  the  plane  (j  =  TJ)  of  maximum  H.  Is in  (j  +  u)   as  a 
reference  plane.   In  Figure  10  below  search  depth 
HTJ  | sin  (TJ  +  u)  |  =  H(TJ+1j|sin  (TJ+1+u)  |  . 

In  this  case  where  more  than  one  of  the  planes  j  has  the  same  maximum 
H.  I  sin  (j  +  p.)  ,  the  reference  plane  is  chosen  to  be  that  with  the 
least  H.  cos  (j  +  \i) .  Above  search  depth,  each  plane  j  such  that 
H.  sin  (i  +  u)<Dr,„  -  CDA  and  i>TL  defines  a  level  L, .  Also  each 

plane  j  such  that  H.  sin  (i+  u)2:  D  -  CDA  defines  a  level  L, » 
J  jXJr/WS  d 


^9 


Referring  to  Figure  11.  the  volume  elements  LADV ,  .  MADV ,  «  and 

du'     du' 

RADV,   (d  =  1,  *";  D;  u  =  2,  •••,  U)  are  significant .  The  subscript  d 

indicates  the  level  L,  and  the  subscript  u  indicates  the  search  turn 

number.  The  level  L-  is  the  final  level  above  the  critical  depth  CDA 

and  L=  =  TES .  The  corresponding  volumes  below  the  search  depth  are 

LBDV  ,   MBDV  ,  and  RBDV   (d  =  1,  •'•,  Dj  u  =  2,  •••,  U  ).   The  scanned 

volume  in  level  L  above  the  search  depth  on  the  u  th  search  turn  is 

ASV .   =  LADV .   +  MADVj   +  RADV,, 
du       du       du       du . 

One  or  more  of  these  volume  increments  may  be  zero  as  is  LADV,-, 

d3 

in  Figure  11.  The  corresponding  scanned  volume  below  the  search  depth 

is  BSV,  . 
du 

Referring  to  Figure  12, 

the  volumes  AMV   ;  d  =  D+l;  '",  d'j  u  =  1,  "',  U  are  significant . 

The  corresponding  volumes  below  the  search  depth  are  denoted 

BUV^  (d  -  D  +  1,  -.,  D1;  u  «  1,  ••.,  U  ) 

For  levels  L,  below  depth  CDA  and  above  search  depth  and  for  the  u  th 
d 

search  turn,  AMV   is  the  trapped  volume  of  target  uncertainty  remain- 

ing  after  the  ping  cycle  is  complete.  The  volume  scanned  in  level  L, 

above  the  search  depth  on  the  u  th  search  turn  is  thus 

ASV_,   =  AMVw   .  N  -  AMVJ 
du      d^u-1;      du 

The  corresponding  scanned  volume  below  the  search  depth  is  BSV,  „ 
In  this  step  all  calculations  are  done  for  ping  zero  and  the 
quantities  computed  for  ping  zero  are  saved  to  be  modified  later  as 
appropriate  for  the  i  th  ping.  The  total  scanned  volume  for  ping 
zero  of  the  u  th  search  turn  is 

50 


TD 

4) 

C 

C 

«8 

ty 

o 

J3 

CO 

U 

>* 

o 

,-t 

u 

T-t 

3 

i 

4J 

o 

u 

TJ 

0 

(0 

4J 

00 

M-l 

cU 

c 

4) 

V 

•H 

■H 

M 

a 

<fi 

o 

4) 

r-4 

J3 

a 

4J 

<! 

CO 

o 

Q 

0 

4J 

o 

e 

i-t 

B 

J3 

f-« 

5 

4J 

o 

a 

« 

T) 

4) 

> 

T3 

« 

60 

J3 

c 

.-< 

1-1 

01 

-o 

a 

O 

4) 

o 

i-l 

e 

«-i 

u 

5 

CO 

•r-l 

as 

M 

U 

CO 

u 

(0 

4) 

d 

4) 

(0 

■H 

> 

OS 

.-l 

o 

4) 

.fi 

Pi 

4) 

a) 

OS 

> 
OS 

T3 

4J 

X. 

h) 

a 

4) 

4) 

f-l 

f-l 

a 

V 

OS 

os 

> 

> 

o 

4> 

•H 

CO 

i~l 

3 

t3° 

4) 

0) 

W 

•U 

r-l 

4) 

a 

60 

0 

« 

M 

0} 

>* 

oS 

0) 

i-H 

4-) 

u-t 

OS 

<4-l 

o 

o 

e 

4J 

CO 

0 

M 

oS 

•H 

4) 

4) 

4J 

> 

M 

U 

< 

a> 

-o 

m 

4) 

0 

a 

• 

CO 

i-i 

4J 

(0 

■-« 

J3 

o 

00 

U 

CO 

i-t 

O 

■H 

M 

4> 

3, 


51 


'I 


CO 

4J 

c 

4) 

CO 

0 

V 

4 

CVJ 

M 

4) 

T3 

a 

u 
< 

i 

4) 
H 

a 

0 

eg 

o 

< 

6 

a 

3 

4) 

o 

i-4 

e 

0 

i-j 

J3 

> 

pH! 

4-» 

Ou 

co 

CM 

0) 

4J 

O 

•o 

c 

0) 

u 

i-H 

03 

J5 

0) 

4> 

60 

U 

)-i 

fl 

•H 

a 

H 

4J 

V 

i-l 

Vj 

O 

)-l 

*J 

O 

o 

w 

« 

3 

ce 

4) 

5 

<: 

l-i 

t-i 

< 

«» 

« 

J3 

4) 

n 

° 

•o 

<! 

ro 

h4 

T3 

r-t 

0 

s 

4> 

1-1 

^ 

> 

•o 

4) 

^» 

4) 

1-1 

§ 

1 

<a 

i— 1 

t-t 

t) 

0 

• 

a 

B 

> 

cu 

6 

3 

s 

« 

i-i 

CO 

3 

CO 

O 

u 

i~l 

> 

e 

o 

U-l 

4) 

> 

o 

CO 

CO 

4J 

4) 

*o 

c 

e 

H 

4> 

o 

4) 

a 

C 

-H 

CO 

4) 

d 

4J 

at 

M 

m 

u 

u 

u 

V 

a 

O 

CO 

09 

u 

fa 

0 

M 

S 

i-t 

(0 

86 

(0 

H 

(0 

3 

o 

V 

•U 

H 

o 

H 

U 

o 

<: 

«< 

c« 

CVJ 

1-1 

4) 

a 

fa 


52 


D'  D" 

SV0u  "    S     ASVdu     ♦         ^      BVdu 
d=l  d=l 

The  output  of  this  step  is  the  number  of  pings  (M)  in  a  search 
turn,  the  maximum  available  number  of  search  turns  U^  several  tables 
of  parameters  described  in  appendix  Ck   used  in  computing  the  special 
volumes,  and  tables  of  the  following  special  volumes  for  ping  zero; 

LADVdu,  MADVdu,  RADVdu;  d  =  1,  ••.,  D;  u  =  2,  •••,  U 

LBDVdu,    MBDVdu,    UWfcj    d   =   1,    •••,    D;    u  -  2,    -..,    U 

AMV      j    d  =  D  +    1,    •  •  •,    D1  ;    u  =   1,    •  <  •,    U 

BMV      ;    d  =  D  +   1,    •  •  • ,    D1  ;    u  =   1,    •  •  • ,    U 

ASVdu;    d  =   1,    •..,    D* ;   u  =  2,    ...,   U 

BSVdu;    d  =   1,    •••,    D!;   u  =  2,    ••-,  0 

SV^;   u  =  2,    ••',    U 
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5.  First  turn  scanned  volume:  i  th  ping 

Throughout  the  search  phase,  the  weapon  is  displaced  from  the 
center  of  the  target  volume  cylinder  a  horizontal  distance  R  .  The  base 
radius  of  the  target  cylinder  is  some  distance  R  at  the  time  the 
weapon  splashes  into  the  water.  The  radius  then  increases  at  a  rate 
(S   )  equal  to  the  target  cruise  speed  to  a  distance  R..  at  the  time 
of  the  first  weapon  ping.  Thereafter  the  radius  increases  at  a  rate 
(S   )  equal  to  the  alerted  target  escape  speed  and  has  a  value  R„ 

id  1 

upon  commencement  of  the  i  th  ping  cycle. 

In  the  situation  of  Figure  13,  a  top  view  of  the  target  cylinder 
of  increasing  volume  with  the  first  few  ping  orientations  and  scanned 
volumes  in  the  first  search  turn,  the  insonified  volume  of  ping  1 
includes  a  large  volume  outside  the  target  cylinder.   Such  volume  does 
not  contribute  to  acquisition  probability  because  there  is  an  assumed 
zero  probability  that  the  target  is  outside.  Therefore  the  total 
insonified  volume  must  be  reduced  to  the  useful  scanned  volume  inside 
the  target  cylinder.   In  the  situation  depicted,  the  amount  of  dis° 
carded  volume  is  reduced  with  each  ping  until  ping  9  when  all  of  the 
insonified  volume  is  within  the  target  volume  cylinder. 

The  projection  of  the  radial  (R.,  ,  p\,  )  in  plane  j  to  the  hori- 

jR    JK 

zontal  plane  through  the  transducer  is  the  line  OP  of  Figure  Ik 
denoted 

R.^    where 

Rjk  =  V  cos  <J  ^  > 


Figure  13:  Top  view  of  target  cylinder  of  increasing 

volume  with  the  insonified  volumes  of  several 
sample  pings  of  the  first  search  turn  depicted 
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Figure  Ik :  The  projection  of  plane  j  in  the  horizontal  at 

search  depth.   Line  00*  is  the  plane  j  reference 

axis  projection.  Line  IP  is  CR.,  .   Line  OP  is 

R1.,  .   Line  IP'  is  R..  Line  OP1  is  r'.  „,  . 
jk  i  ijk 
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Use  of  the  law  of  cosines  enables  determination  of  the  distance 
I  P  denoted  CR.,  .  This  distance  is  compared  to  the  radius  (R. )  of  the 
target  cylinder  on  the  i  th  ping.   If  the  radius  (R,  )  is  smaller^  it 

X 

defines  the  new  triangle  IOP  and  the  projection  of  the  side  OP  back  to  tl 
plane  j  defines  a  new  acquisition  range  (R, ..)  to  be  associated  with  the 

X  J  K» 


angle  B.,  .   That  is 

ijk 
ijk  "   cos  ( j  +  u  ) 


However  if  the  distance  CR.,  is  the  lesser  of  the  two  then 

R .  .,  =  R , , 
ijk    jk 

The  remaining  portion  of  the  computation  for  scanned  volume  (SV. .) 
on  the  i  th  ping  of  the  first  search  turn  is  identical  to  that  in 
step  3  for  SV  ,. 

The  output  of  this  step  is  the  table  of  volumes  SV.,j  i  ^  1,  °  °  °  9   M= 

and  the  distances  R  and  R.:  i  =  1,  °°'*   M- 1 . 

o      l       *     ' 
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6.   Scanned  volumes  subsequent  to  first  turns   i  th  ping 

Just  as  the  first  turn  scanned  volume  had  to  be  corrected  for  the  i 
ping  because  of  the  effect  of  the  target  cylinder  radius,  so  must  the 
scanned  volumes  on  subsequent  turns. 

Define  the  intersection  of  the  ping  i  bisector  with  the  target 
cylinder  radius  at  transmission  of  ping  i  to  be  a  distance  CH  from  the 
weapon. 

Referring  to  the  vertical  cross=section  of  the  insonified  volume 
in  Figure  15,  verticals  through  four  of  the  possible  intersections  have 
been  plotted  and  labeled  CHI  through  CH4  „  The  effect  of  each  of  these 
four  cases  will  be  discussed  as  it  applies  to  a  level  above  and  a  level 
below  the  critical  depth  CDA  (i.e:  levels  2  and  U).  The  treatment  of 
the  situation  below  the  search  depth  is  identical . 

In  the  case  of  CHl  the  first  turn  scanned  volume  is  reduced., 
However  the  scanned  volumes  for  all  subsequent  turns  are  the  same  as 
those  computed  in  step  four  for  the  general  ping.  This  is  because 
CHl  lies  to  the  right  of  all  areas  representing  equivalent  volumes . 

In  the  case  of  CH2  there  is  still  no  reduction  for  volumes  in 
level  2.  However  there  is  a  small  reduction  in  the  volume  scanned  on 
the  second  search  turn  in  level  h .  The  volume  reduction  is  that 
volume  associated  with  the  area  in  level  k   to  the  right  of  CH3» 

It  is  computed  as  a  reduction  in  AMV,  .  (the  volume  missed  in 
level  Li  above  search  depth  on  the  first  search  turn). 

Thus  (ASV22)ping  .  equals  (ASV^,)^  Q  but 

<ASVplng  1   1S   USS   than  (AS\2>Plng  o- 
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In  the  case  of  CH3  there  is  a  reduction  in  the  scanned  volume  on 

the  second  search  turn  in  both  levels  L^  and  L,  „   In  level  L_  the  volume 

2      4  2 

(ASV\^)  .    .  is  less  than  (ASV^^)  .    .  by  an  amount  equal  to 
v    22 'ping  l  N    22 'ping  0  ' 

RADV   plus  a  portion  of  MADV   . 

In  level  L.  the  reduction  still  amounts  to  a  portion  of  AMV.  . . 

In  the  case  of  CH4  there  is  no  scanned  volume  in  level  L_  above 

2 

search  depth  on  the  second  turn.   That  is  ASV   =  0  for  ping  i „ 
Further  there  is  a  small  reduction  in  level  L_  in  the  scanned  volume 
for  the  third  search  turn.   For  level  L,  there  remains  some  scanned 
volume  on  the  second  turn  and  there  is  also  a  reduction  in  scanned 
volume  for  the  third  turn. 

The  output  of  this  step  is  a  table  of  scanned  volumes 
SViu;  i  =  1,   ••,  M-l;  u  =  2,  ••-,  U. 
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7.  Target  Volume:  i  th  ping 

Consider  that  portion  of  the  target  cylinder  volume  computed  at 

a  base  radius  of  R.  within  the  radials  bounding  the  ith  ping.   See 

Figure  16.   The  area  contained  between  the  boundary  radials  of  the 

i  th  ping  and  within  the  cylinder  of  radius  R.  multiplied  by  the 

height  of  the  cylinder  is  this  volume.  The  probability  that  the 

target  is  in  this  i  th  volume  slice  of  the  target  cylinder  is  the 

volume  of  the  slice  divided  by  the  total  cylinder  volume „   This 

probability  (PQ.;  i  =  1,  •••,  M-l)  remains  constant  throughout  the 

problem  because  of  the  opening  target  assumption  number  six0 

From  Figure  16: 

area  of  horizontally  lined  region 
1  TT  R  2 

The  target  volume  (TV.j  i  =  1,  °,of   M~l)  associated  with  the  i  th 
ping  on  the  first  search  turn  lies  between  the  same  radial  limits  but 
within  a  target  cylinder  of  radius  R.  instead  of  R. .  That  is 

TV.  =  (area  of  vertically  lined  region)  D 

Figure  16  shows  the  angle  ( b'  )   assigned  to  the  M  th  ping  to  be  less 
thanb  .   It  is 

b  '  =  360  -  (M-l)b 

This  smaller  angle  is  required  to  complete  the  search  turn  so 
that  the  i  th  ping  of  the  u  th  turn  is  oriented  in  the  same  direction 
as  the  i  th  ping  in  the  (u  +  1)  th  turn.  Thus  the  M  th  ping  is  in 
part  ficticious. 

The  output  of  this  step  is  thus  the  tables  of  TV.;  i  -  1,  ° » ° f   M~l 
and  of  PQ.;  i  =  1,  •••,  M-l  and  the  angle  b"  . 
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ping  1 


Figure  16:  Top  view  of  target  cylinder  of  increasing  volume  „ 

Lines  CT  and  CT1  are  boundary  radials  of  ping  i„ 

p 
Horizontally  lined  area  is  PQ.   times  "^R^ 

Vertically  lined  area  times  cylinder  height  is  TV 


I 
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8.  Acquisition  Probabilities:  I  pings 

The  instantaneous  probability  of  acquisition  on  the  i  th  ping  is 
simply  the  ratio  of  the  volume  effectively  scanned  to  the  volume  of 
target  uncertainty. 
Define:   g.   =  P  ["acquisition  during  the  i  th  ping  cycle  of  the  u  th 


search  turnj 


Then     glu  =  SViu  P^  /  TT±  J  i  =  1,  •",  M-l$  u  =  1,  •••,  U 


Recall  that  PQ.  is  the  probability  that  the  target  is  within  the 
radial  limits  of  ping  i. 

The  volume  TV.  is  the  total  volume  of  target  uncertainty  within 
the  radial  limits  of  ping  i.  The  volume  SV„.  is  the  volume  scanned 
on  the  first  turn.  The  volumes  scanned  on  subsequent  turns  are  volumes 
measured  inside  the  first  turn  scanned  volume  and  under  the  conditions 
existing  for  the  tirst  turn.  These  volumes  are  then  moved  intact  with 
the  opening  taget.  Thus  the  applicability  to  all  turns  of  the  volume 
TV.  measured  using  the  first  turn  cylinder  radius „ 

The  instantaneous  probability  of  detection  on  the  M  th  ping  is 
approximated  by 

«*,  =  «<*.«»  &*>*>♦  «i<u*i)  »-*>/& 

The  M  th  ping  is  alloted  the  same  time  interval  as  any  other  ping0 
Errors  involved  should  be  small. 
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There  remains  the  requirement  to  relabel  the  g,   to  be 

iu 

g.;  i  ss  1>  ' "  f    I>  to  compute  I  as  a  function  of  w.  as  described  in 
the  beginning  of  the  chapter,  to  compute  P  as  a  function  of  I,  and 

A 

to  tabulate  the  values  of  P  versus  w.  as  w.  is  allowed  to  vary  from 
W  to  W  .  This  final  table  is  the  input  to  the  dynamic  programming 
model . 
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9.      Input  Requirements 

A.  Medium 

c  =  velocity  of  sound  in  water  (yd/ sec) 

a  =  attenuation  coefficient  (db/yd) 

10  log  m  =  surface  reverberation  scattering  coefficient  (dl 
s 

10  log  m^  =  volume  reverberation  scattering  coefficient  (db 

B .  Target 

TS  =  target  strength  (db) 

TN  =  radiated  noise  level  of  the  target  (db  at  1  yd) 

D_.  =  maximum  target  depth  (ft) 

TM 

S   =  target  speed  during  weapon  search  (yd/sec) 
S   =  target  cruise  speed  (prior  to  becoming  alerted  on 
the  first  ping  cycle)  (yd/sec) 

C .  Weapon 

WS  =  transmitted  on=axis  source  level  of  weapon  (db  at  1  y< 
AT  =  ping  length  (sec) 
PRR  =  weapon  active  sonar  pulse  repetition  rate  for 

search  (pings/sec) 
Af  =  receiving  bandwidth  (cycles/sec) 

WEL,  =  weapon  enabling  level  tor  a  range  of  lOuk  yds  (dbj 
SN  =  equivalent  isotropic  self~noise  level  at  search 

speed,  deptn,  and  turn  rate  (dbj 
RD  =  recognition  differential  for  5@$  probability  of 
acquisition  (db) 
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b.   =  normalized  pattern  function  in  the  vertical  plane 

for  transmitting  pattern  (a  numeric) 
b   =  normalized  pattern  function  in  the  vertical  plane  for 

receiving  pattern  (a  numeric) 
b   =  normalized  pattern  function  in  the  azimuth  plane  for 

transmitting  pattern  (a  numeric) 
b,   =  normalized  pattern  function  in  the  azimuth  plane  for 

receiving  pattern  (a  numeric) 
J   =  limiting  angle  for  off<=axis  reception  in  the  vertical 

plane  (deg) 
«,  <*•  =  limiting  angles  for  off~axis  reception  in  the  azimuth 

plane  (deg) 
DR  =  directivity  index  for  receiving  (db) 
u   =  angle  of  attack  (pitch)  of  weapon  (deg) 
0)_  =  search  rate  of  turn  (deg/sec) 

0 

DS  =  angle  of  dive  in  diving  to  search  depth  (deg) 

D„.  =  weapon  search  depth  (ft) 
WS 

S   =  weapon  speed  in  dive  to  search  depth  (yd/sec) 

S  „  =  weapon  speed  during  search 
WS 
I 

W   =  maximum  allowable  weapon  fuel  weight  (lb) 

w.   =  weapon  fuel  expenditure  rate  at  search  depth  and 

speed  (lb/sec) 
w   =  average  weapon  fuel  expenditure  rate  in  dive  to 

search  depth  (lb/sec) 
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D.  Delivery  weapons  system  capability 

B  =  radius  of  target  volume  cylinder  at  time  of  weapon  splash 
Statistical  averages  and  distributions  are  available  for  the  medium 
inputs.   Intelligence  reports  and  estimates  or  CNO  requirements  may 
provide  the  target  inputs  except  for  the  target  speeds  which  may  be 
selected  in  the  higher  ranges  in  accordance  with  criterion  one. 
The  weapon  inputs  are  available  from  weapon  designers  at  least  in  the 
form  of  educated  estimates.   Finally  the  input  R  may  be  estimated 
from  observations  of  the  performance  of  available  delivery  weapons 
systems . 
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CHAPTER  IV 
PROBABILITY  OF  HIT 
Consider  a  tactical  situation  in  which  a  weapon  acquires  a  submarine 
target  at  some  initial  slant  range  and  make  the  following  assumptions. 
Assumption  1:  The  weapon  makes  an  attack  and  a  series  of  independent 
reattacks  on  the  target  culminating  in  either  a  hit  or  fuel  exhaustion. 
Assumption  2:  At  the  time  of  acquisition,  the  target  is  opening  at 
maximum  speed.  Target  course  and  speed  are  constant  throughout  the 
attack  and  reattack. 

Assumption  3:  On  the  attack  and  after  reacquisition  on  each  reattack 
the  weapon  maintains  contact  and  closes  at  attack  (maximum)  speed. 
Assumption  k-i     The  target  operates  in  depth  within  specified  minimum 
and  maximum  limits  such  that  the  maximum  depth  limit  is  not  greater 
than  the  weapon  maximum  operating  depth. 

Assumption  5:  The  depth  at  which  a  hit  may  take  place  is  a  random 
variable  denoted  the  hit  depth  which  has  an  assumed  probability 
density  function. 

Assumption  6:  The  average  weapon  depth  for  the  attack  is  the  arithmetic 
average  of  the  search  depth  and  the  average  hit  depth.  The  weapon  and 
target  depths  for  the  final  stage  of  the  attack  and  throughout  each 
reattack  are  equal  at  some  depth  within  the  limits  of  assumption  four. 
Assumption  7:  The  attack  and  each  reattack  is  from  astern  the  target. 
There  is  some  close-in  range  at  which  the  weapon  sensors  lose  contact 
with  the  target.  On  losing  contact  with  the  target  the  weapon  continues 
on  a  straight  track  at  attack  speed  for  a  preset  dead  time.  Upon  com- 
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pletion  of  the  dead  time  run,  the  weapon  instantaneously  slows  to 

search  speed  and  commences  a  left  turn  at  reattack  turn  rate  and 

constant  depth. 

Assumption8:  The  weapon  reattack  turn  is  a  circle  with  no  advance  or 

transfer.  The  weapon  reacquires  the  target  upon  completion  of  one 

complete  reattack  turn  and  instantaneously  increases  to  attack  speedy 

Assumption  9 J   At  the  time  the  close-in  weapon  sensors  lose  contact 

with  the  target,  there  exists  an  intersection  of  the  weapon  track 

extension  with  the  vertical  plane  through  the  submarine  target  center 

which  is  perpendicular  to  any  plane  containing  both  the  target  and  the 

weapon.   Let  the  target  projection  in  the  vertical  plane  be  approximated 

by  a  circle  denoted  the  target  circle  centered  on  the  target  longitudinal 

axis  with  a  diameter  equal  to  the  length  of  the  beam  dimension  of  the 

target.   Let  the  circle  center  be  the  origin  of  a  doubly  infinite 

axis  containing  the  intersection.  Assume  that  the  intersection  co= 

ordinate  on  the  axis  is  a  normally  distributed  random  variable  with 

mean  zero  and  standard  deviation  a    . 

Assumption  10:  A  hit  occurs  if  the  weapon  track  at  the  closest  point 

of  approach  to  the  target  center  passes  within  a  distance  of  the  center 

not  greater  than  the  radius  of  the  target  circle  plus  a  small  proximity 

correction  factor. 

Define: 

P     =  P     (N)  =  P  [hit   in  N  attacks] 

h     =  P  ["hit   on  the  n  th  attack]    j  m  *  1,    •••,    N 


PH  (0)  =  0 
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Then 

N 

p  =  i  -  TT    (i  -  h  ) 

H      n=l       n 
The  derivation  is  analogous  to  that  for  acquisition  probability 0 

Define: 

t  =  time  the  weapon  sensors  lose  contact  with  the  target  on  an 

attack  or  reattack  (zero  reference  time) 

RH  =  range  from  the  weapon  to  the  target  center  at  time  t  (yd) 

t    =  time  from  t   of  closest  point  of  approach  of  weapon  and 

target  (sec) 

MRH  =  distance  traveled  by  weapon  from  time  t   to  time  t    (yd) 

O  \w  c£\ 

S   =  weapon  attack  speed  (yd/sec) 

WA 

S_.  =  target  maximum  speed  (yd/sec) 
TM 

BW  =  target  beam  width  (yd) 

VT  =  weapon  proximity  correction  factor  (yd) 

r   =  coordinate  of  intersection  of  the  weapon  track  extension 

with  the  vertical  plane  (A)  through  the  submarine  target  center 

perpendicular  to  any  plane  containing  both  target  and  weapon  at 

time  t  (yd) 
o  w  ' 

r'   =  coordinate  of  intersection  of  the  weapon  with  the  vertical  plane 
(A*  )  parallel  to  plane  A  containing  both  the  submarine  target 
center  and  weapon  at  time  t    (yd) 

\s  in 

F  ix)  =  cumulative  normal  probability  distribution  with  mean  zero  and 
standard  deviation  one  evaluated  at  the  point  x 
Consider  the  attack  situation  depicted  in  Figure  17 . 
The  xy  plane  is  the  vertical  containing  the  weapon  and  target  at  time 
t  .   Planes  A  and  A1  are  parallel  to  the  yz  plane  and  contain  the 
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bases  of  two  right  circular  cones  with  a  common  apex  at  the  origin „  The 
line  OR  is  the  weapon  track  and  its  orientation  with  respect  to  the  x 
axis  determines  the  orientation  of  the  axes  LL  and  L  L  on  which 
coordinates  r  and  r  are  measured.  The  common  intersection  of  the  cone 
of  base  radius  r' ,  the  plane  A* ,   and  the  xy  plane  is  a  line  of  length 

2|r'|  .   By  similar  triangles 

i     MRH 
r  =  r 

RH 

Since 


MRH  =  KH  SWA  /  (SWA  -  STM) 


then 


r'  =  r  SWA  /  (SWA  "  STM> 


I 


Thus  r   is  a  random  variable  with  mean  zero  and  standard  deviation  0 


a'  =  a  SWA  /  (SWA  -  STM) 


Therefore 


h     =  P 


[2|r.| 


BW  +  2  VT 


] 


=  P  [|r«|  <(BW/2)   +  VT] 


=  2  F 


(BW/2)   +  VT 


<*> 
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b  =  2  F  . 

n     0 


(BW/2)  +  VT 


*  SWA  /  <SWA 


stm> 


Define; 

HD  =  hit  depth  (ft) 

f  (x)  =  probability  density  function  of  the  random 
HD 

variable  hit  depth 

D_..  =  maximum  target  depth  (ft) 
TM 

MIN  =  minimum  target  depth  to  be  considered  (ft) 

Define  a  set  of  Q  equal  depth  subintervals  in  the  depth  interval 

D   ,  D   J  such  that  the  depth  at  the  lower  boundary  of  the  q  th 

11  isDMiN  +  q  pWSmi]- 

P   =  P  hit  depth  is  in  the  q  th  depth  sub interval! ;  q=l,  ° • ° ,  Q 
Assume  as  in  Figure  18: 


subinterva! 


Define 


f-nW  =  2  (x  -  Dm)  /  (D 


■HD 


MIN' 


TM   WMIN 


)' 


Then: 


Eq 


_  2q  '  1 


•  q  =  I   ...,  Q 


Define: 

i 


w   =  fuel  used  in  the  n  th  attack  or  reattack  (lb) 


2n 


=  initial  acquisition  range  (yd) 
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Figure  18:   Assumed  probability  density  function  of  the 
random  variable  hit  depth. 


Ik 


w   =  weapon  fuel  expenditure  rate  at  attack  speed  and 

at  the  average  depth  for  the  attack „   (lb/sec) 
HD  =  average  hit  depth  (ft) 
Then 

55  "  2/3  DTM  +  l/3  °M1N 
and  the  average  depth  for  the  attack  is 

Vll!  =  1/2  °ws  ♦  1/3  DTM  +  1/6  DWN 

For  the  attack 

R   •  w 


21    SWA  -  STM 


Define: 


w    =  weapon  fuel  expenditure  rate  at  search  speed  and  at 
the  midpoint  of  the  q  th  depth  interval,  (lb/sec) 

w    =  weapon  fuel  expenditure  rate  at  attack  speed  and 

at  the  midpoint  of  the  q  th  depth  interval,  (lb/sec) 

t„  =  dead  time  from  t  at  which  the  weapon  commences  the 
DT  o  r 

reattack  turn  (sec) 
t    =  time  from  t  at  which  the  weapon  completes  the  reat 

turn  and  reacquires  the  target  (sec) 
Ol.   =  reattack  rate  of  turn  (deg/sec) 


Then 

• 

0 

tCT   = 

trvm     + 

DT 

360° 

tCPA 

RH 

SWA  " 

STM 
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The  weapon  will  spend  the  following  time  at  attack  speed  during  each 
reattack:   See  Figure  19. 


ta  =:  St  "  Spa  + 


^TCT  "  SpA^  STM°  ^St'SpA'  SWA 
SWA  '  STM 


The  weapon  will  spend  the  following  time  at  search  speed  during 
each  reattack: 


Ts  = 

360° 

The- 

refore: 

Q 

<       ■ 
2n 

■>  I 

q=l 

n  =  2 

,     •••,    N 

T.  w_   +  T_  vnj_ 
A  Atq    S   Stq 


Eq 


Define : 

S   ■  minimum  preset  value  of  attack  speed  (yd/sec) 
o 

w'   =  available  fuel  weight  for  hit,  event  2,  (lb) 

w' '  =  available  weight  for  modifications  to  the  propulsion 

system  to  provide  a  maximum  speed  capability  in  excess 

of  S  (lb) 
o 

N  is  a  function  of  w'  according  to 

N+l 


^>~  w'  <w'  <^>   wL 

/  2n    2  ^  y        2n 


n=0 


n=0 


However,  recall  that  w_  is  a  function  of  attack  speed  which  is  it= 

'  2n  r 

self  a  function  of  w' ' . 
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TPA 


'CPA 


^St^cfa^  ^3M 


(V^i^  s 


DT  "CPA'     WA 


fteT"tePA)  stm 


■> 


^DT^OFA^  SWA  *  (tGT"tCPA)  STM"(tDT"tCPA)SWA  SWAffiwA~STM} 


Figure  19:   Geographic  plot  of  weapon  and  target  tri 

separated  to  show  time  differences.  The 

distancesbelow  the  target  track  are  measured 

from  the  joint  position  at  time  t    . 

UJfA 
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Thus 

N(w»,  v")      ^ 


Pff(N)  =  1  -  (T 

n=l 


1  =■  h    Cw9 ' 
J.       nn^w2 


Recall  w     =  w*    +  w1 ' 
2         2         2 

i       i 

W,  W  =  lower  and  upper  limits  respectively  on  the  range  of  fuel 

weight  (lb) 
W  ,  W   =  lower  and  upper  limits  respectively  on  the  range  of 
weight  for  modifications  to  the  propulsion  system  to 
provide  a  maximum  speed  capability  in  excess  of  S  (lb) 
For  each  value  in  the  range  for  w_ 
min  -f  W1,  W"i<w  <min  -[  W  +  W  ,  Wj* 

the  function  P  is  evaluated  for  each  possible  combination 
H 

of  the  variables  w'  and  w''  as  they  are  allowed  to  vary  such  that 

they  remain  within  their  respective  ranges 
i     i    _i 

W  <wr  <W 
—     2 

ii    ii   _i i 

W  <w  2~W 

and  their  sum  is  equal  to  w  .   Let  K  be  the  number  of  values  of 
w?.  This  procedure  provides  J,  sets  of  values 

W2jk'  W2jk'  PHjk'  j  =  l>    °°°>    Jk 
for  the  k  th  value  of  w  ,  .   (k  =  1,  • • ■ ,   K) 
The  entire  set  of  sets  of  values 

W2k'  V2jk'  W2jk'  PHjk;  j  "  l'    *'  Jk 
k  =  1,  '••,  K 

78 


is  the  input  to  the  dynamic  programming  allocation  model . 

Input  Requirements 

1.  Target 

BW  =  target  beam  width  (yd) 

D   =  target  maximum  depth  (ft) 

S   =  target  maximum  speed  (yd/sec) 

2 .  Weapon 

S   =  weapon  search  speed  (yd/sec) 

W  o 

VT  =  weapon  proximity  correction  factor  (yd) 

CBu  =  weapon  reattack  turn  rate  (deg/sec) 

i  _i 

W,  W  =  lower  and  upper  limits  respectively  on  the  allowable 

range  of  fuel  weight  (lb) 

S   =  minimum  allowable  weapon  attack  speed  (yd/sec) 

II   _i  i 

W  ,  W   =  lower  and  upper  limits  respectively  on  the  range 

of  weight  for  modifications  to  the  propulsion 
system  to  provide  a  maximum  speed  capability  in 

excess  of  S  (lb) 

o  N   ' 

G,  =  Relationship  between  increased  weapon  speed  above  S  and 
increased  weapon  weight 

G_  =  Relationship  between  fuel  expenditure  rate  and  speed  and 
depth.  Recall  that  this  relationship  may  vary  as  a 
function  of  attack  speed  only  in  the  range  for  speeds 
above  search  speed.  This  relationship  is  required  to 
specify  w   ,  w   ,  and  wgt  (lb/sec) 
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3.  Miscellaneous 

RH  =  range  from  the  weapon  to  the  target  center  at  the  instant 
(t  )  during  close-in  attack  when  the  weapon  sensors  lose 
contact  with  the  target  (yd) 
O  =  standard  deviation  of  the  normally  distributed  random 

variable:  coordinate  of  intersection  of  the  weapon  track 
extension  with  the  vertical  plane  (A)  through  the 
submarine  target  center  perpendicular  to  any  plane  con= 
taining  both  target  and  weapon  at  time  t   •  (yd) 

D    =  minimum  target  depth  to  be  considered  (ft) 

M1N 

R  =  initial  acquisition  range  (yd) 

Q  =  number  of  depth  subintervals  in  the  depth  interval 

W  dtm] 

Intelligence  reports  and  estimates  or  CNO  requirements  may 
provide  the  target  inputs.  The  weapon  inputs  are  available  from 
weapon  designers  at  least  in  the  form  of  educated  estimates.  The 
inputs  RG  and  0   are  functions  of  the  medium^,  the  weapon  sonar  and 
guidance  systems,  the  target  physical  dimensions,  and  the  orientation 
relationships  between  the  target  and  weapon.   In  the  absence  of  a 
specific  theory  these  must  be  handled  by  simulation  or  sensitivity 
analysis.  The  input  D is  a  function  of  the  capability  of  the 
weapon  delivery  system  and  the  tactical  situation  as  seen  by  the 
target  commander.   In  the  absence  of  specific  considerations  cri- 
terion one  would  determine  D-,_xl  to  be  zero.   The  input  R.  should  be 

MIN  A 

8o 


selected  close  to  the  maximum  weapon  capability  in  accordance  with 

criterion  one.  The  function  P  may  be  very  sensitive  to  changes  in 

n 

R  and  an  examination  of  results  over  a  range  of  values  for  R  may  be 

A  A 

appropriate.  The  input  Q  should  be  large  enough  to  point  up  the 
effect  of  an  irregularly  shaped  fuel  expenditure  rate  curve  but  small 
enough  not  to  excessively  prolong  the  computer  program „ 
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CHAPTER  V 

PROBABILITY  OF  KILL  DAMAGE 
Consider  a  tactical  situation  in  which  a  weapon  explodes  within 

a  specified  proximity  of  some  portion  of  a  submarine  target  and  make 
the  following  assumptions. 

Assumption  1:   The  hit  depth  is  a  random  variable  which  has  an 
assumed  probability  density  function. 

Assumption  2:   The  point  on  the  target  at  which  the  hit  occurs  de- 
noted the  hit  point  is  a  random  variable  which  has  a  uniform  distri- 
bution  over  the  entire  surface  area  of  the  target ,  The  weapon 
approaches  the  target  from  astern,  loses  contact  at  some  close-in 
range  and  continues  on  a  straight  track.   Since  the  target  may  be 
in  a  turn,  climb  or  dive,  all  points  of  the  target  are  available 
as  hit  points. 

Assumption  3:   The  target  operates  in  depth  within  specified  minimum 
and  maximum  limits  such  that  the  maximum  depth  limit  is  not  greater 
than  the  weapon  maximum  operating  depth.   Divide  the  target  surface 
area  into  P  area  elements,  each  of  which  has  a  constant  hull 
vulnerability  throughout.  Divide  the  interval  of  possible  hit 
depths  into  Q  equal  depth  subintervals  such  that  the  depth  at  the 
lower  boundary  of  the  q  th  subinterval  is 


DMIN  +  < 


D  -D 
TM  MIN 
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where 

D_w  =  maximum  target  depth  (ft) 

D    =  minimum  target  depth  to  be  considered  (ft) 

Define:   P  =  P  (w  )  =  P   kill  damage  given  an  explosive  weight  w  1 

f  (x)  =  probability  density  function  of  the  random 
rlD 

variable  hit  depth 
P   =  P   hit  depth  is  in  the  q  th  depth  subintervalj  j 

q  =  1,  •••,  Q 
P   =  P   hit  point  is  in  the  p  th  surface  area  element  j 

p  =  1,  ...,  P 
HV  (w  )  =  hull  vulnerability  for  kill  damage  of  the  p  th 

target  surface  area  element  at  the  midpoint  depth 
of  the  q  th  depth  subinterval  given  an  explosive 
weight  w  .  This  number  is  one  if  kill  damage  occurs 
under  the  specified  conditions  and  is  zero  if  kill 
damage  does  not  occur. 

w  =  explosive  weight  (lb) 

tit   _i i i 

W   ,  W    =  lower  and  upper  limits  respectively  on  the  range 

of  explosive  weight  (lb) 
Assuming  as  in  Chapter  IV 

,  ,  ,      2  (x-"W 


Otm-  dmin>2 
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Then 


_      area  of  p  th  surface  area  elementj  p  <■  1,  •",  P 
Ep       total  target  surface  area 


P      Q 

(W3)   "  2    ]>     HVpq  ^  PEP  PI 
p*l     q=l 


t i g     _l  g  g 
The  table  of  P  versus  w  as  w_  varies  from  W    to  W 

is  the  input  to  the  dynamic  programming  model . 

Input  Requirements 

1.  Target 

P   =  P  Thit  point  is  in  the  p  th  surface  area  elementj  j 

p  =  1,  ..-,  P 

D_w  =  maximum  target  depth  (ft) 
TM 

2.  Weapon 

iii   _i  i  i 

W   ,  W    =  lower  and  upper  limits  respectively  on  the  range 

of  explosive  weight  (lb) 

3.  Target  and  Weapon 

HV  (w., )  =  hull  vulnerability  for  kill  damage  of  the  p  th 
pq  3 

target  area  element  at  the  midpoint  of  the  q  th 
depth  subinterval  given  an  explosive  weight  (w_) 
k.     Miscellaneous 

D    =  minimum  target  depth  to  be  considered  (ft) 

MIN 

Intelligence  or  CNO  requirements  may  provide  D_„,  P   and  some 

TM    Ep 

of  the  information  for  HV  (w  )  .  Weapon  designers  may  estimate 

mi    _i 1 1         pq  j 

W    and  W    and  weapons  engineers  may  provide  extrapolated 


estimates  of  HV  (w_). 

pqv  3' 


Qk 
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APPENDIX  A 
FLOW  CHART  NOTATION 
The  following  notation  is  used  in  the  flow  charts  of  the  remaininj 
appendices . 

Function  Box 


A  -^>B 


The  function  box  contains  a  number  or  expression  to  be  evaluated 
A  which  is  to  be  placed  in  memory  cell  B. 

Choice  Box 


The  choice  box  contains  a  comparison  of  the  number  or  expres- 
sion A  to  the  number  or  expression  B.   Exits  are  chosen  as  follows 
Exit  Label 


< 
< 
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In  Box 

INPUT 


The  in  box  contains  an  input  to  the  program .   In  appendices 
C  1  through  C  8  the  output  from  one  appendix  is  carried  over  intact 
into  the  next  appendix.   Inputs  to  be  used  for  the  first  time  in  the 
sequence  of  flow  charts  are  placed  in  a  separate  in  box. 

Out  Box 

OUTPUT 


The  out  box,  one  for  each  appendix,  contains  the  output  of  the 
appendix  program. 

Index  Box 


1  +  i 


The  index  box  contains  a  manipulation  on  an  integer  in  an 
index  register  and  is  the  tool  which  provides  for  loops  in  the 
program. 
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Subroutine  Box 


(RJk 


The  subroutine  box  indicates  a  return  jump  to  the  subroutine 
numbered  by  k.  The  flow  chart  of  each  subroutine  is  shown  after 
the  flow  chart  of  the  first  routine  using  it. 

Return  Box 


RTN 


A  return  box  indicates  the  end  of  a  subroutine  implying  a  jump 
back  to  the  main  routine. 

Fixed  Connector 


d> 


A  fixed  connector  is  used  to  connect  remote  parts  of  the  flow 
chart  with  one  another. 

Assertion  Box 


Note  j 


An  assertion  box  refers  the  reader  to  a  note  found  at  the  end  of 

the  flow  chart. 
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Instruction  flow  is  generally  from  left  to  right  and  from  top  to 
bottom.  All  boxes  are  entered  from  the  left  or  from  the  top  and  all 
exits  from  boxes  are  made  down  or  to  the  right.  The  arrow  is  used  to 
further  indicate  the  sequence  of  computer  instructions.  Arrows  are 
used  at  all  line  intersections  to  indicate  entry  into  the  program 
between  two  boxes . 
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Note  1 


APPENDIX  B 
THE  ALLOCATION  MDDEL 

INPUT 


_L 


PHJk^v2jk»  w2jk^  w2jk»  w2jk°  *  S  lf 


•  o  o        j    * 
9    «1«>J> 


k  ^  1,  •••,  K 


•^(w^t  gj1  (wj»)»  g3(w3)p  Aw,  X 
WnJ  n  «  1,   "%  H 


*N^*MAX 


^i^iV 


Wll^E-*W 


■IE 


-^•D 


L(g) — (i-»k> 


W2k-^W 


5. 


d?-^B 


0_>WJ(W) 


0-»w'«(W) 


-<w> 


n 


In 


V  -  X    ei^V  *  g2*  <V2}k>       ♦  *1    W  -  w2jk  *  w2jk 


a  :  b 


> 


A-^B 


W2jk—  C 


Mil       ^C» 


-Mi  ♦  *  -^i 


fr  *j^j>t 


1  ♦  k-^*k 


Note  5 


W  L= L 


Note  6 


3 


-|B->h2(v) 
C-»wj(W) 
C«— >wJ»(V) 


Note  7 


W  ♦  Aw— =>W 


a^ivA 
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s- 

\ 

WMAX^D 

w31  -^  W  ^  E 

J 

(rjt^ 

/■%  _=>>  «\ — , 

KZ 

Note  8 

Wn      >  WMAX 
fc    (W  )  — >*  F» 

VW-**Pn 

WMM-Pn^^WMM 

y  n' 

n 

u    =  G  <=  H 

MAX    n   nn 


W 


MAX 


P 


"lOW-*-^ 

S3(Pn'^'3n 

v,  +  v'  +  v'  ■  +  v~ 

In    2n    2n    3n 


n 


V*"-»»n)- 


OUTPUT 


Fn»  Gn»  V 
\  V  Pn»  V 


n  -  1, 
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Note  2 


RJ1) u 


«co-e*B 


0^wi(W) 


E->V, 


In  PA(w1)  -  ^(v^)— >A 


In  PE(w3)  -  X«3(w3)  ♦  h2(W  «  w^)— ^A 


A-^B 


»r^c 


w„   ♦  Aw-^-Wo 


facte  3 


B-^h^W) 
C-^W1(W) 


RTN 


Note  U 


i 


0 


<D 


W  ♦  Aw  -^  W 


Note  1:   Provision  is  made  for  a  set  of  payload  weight  limits  (w  )  ii 
case  it  is  desired  to  examine  a  range  of  values.  The  input 


Note  2: 


W„  is  to  be  the  greatest  of  the  W  . 
N  n 

The  subroutine  RJl  computes  the  h.(y)  and  the  w. (w)  for 
i=l,  3«  h  is  handled  separately. 
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Note  3 J  This  loop  provides  that,  for  the  W  considered,  the  function  h„ 

is  maximized  over  the  w. .  For  h.  the  maximum  meaningful  value 

in  both  loops  is  w  =D  so  the  subroutine  is  stopped  by  the 

1R 

comparison  W  :  D.  However  for  h  ,  the  maximum  meaningful 
value  in  the  two  loops  is  different,  the  one  in  this  loop  being 

w__  and  the  one  in  the  other  loop  being  Ww.„  =  D. 

jo  MAX 

Note  k.      This  loop  stores  the  h  (w)  and  w.(w)  obtained  as  W  is  allowed 

to  vary  over  its  range  of  meaningful  values. 
Note  5:  This  loop  maximizes  the  w  .,  over  j  for  a  given  k. 
Note  6:  The  maximum  meaningful  value  of  w  for  a  given  W  is 


min  -Tw2K>  W  I 


»} 


Note  7:  The  h  (w),  w'(w),  and  w'*(w)  are  stored  for  each  value  of  Wu 

Note  8:  This  loop  works  backwards  through  the  stages  to  retrieve 

the  optimum  weight  values  for  each  input  pay load  maximum „ 

Note  9:   The  V  are  the  optimum  volume  sums  to  be  compared  with  the 

input  payload  maximum  allowed  volume  Vo   If  the  solution  is 

not  optimal  with  respect  to  the  restraints,  a  new  X  is  chosen 

and  the  program  repeated.   If  the  final  solution  is  obtained^ 

then  fuel  weight  is  G  +  I  ,  weight  for  added  attack  speed  is 

H  and  explosive  weight  is  P  . 
n  n 
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APPENDIX  CI 
AVAILABLE  ECHO  EXCESS  (fyk) 

INPUT 


c,  a,  10  log  in  j,  TS»  TN9  WS9  AT9  PRR  «  Afs 
Dyg»  J»  SN,  RDP  blV  bgj  b«e  b.P  DR9  ji9  o^ 

i  WEI^J  k»l,  •••,  K« 


fHoto  2 


SN  ♦  10  log  Af  °>  PR  ♦  RD-^>M(SN 


(Not*  3 


Lj(l/2)    (WS  ♦  tS  -  M(SN))-^H 


XO-^k^, 


L^klT 


5> 


Not*  U 


120  log  100k  »  IQOak— >Hfc — [FL    J  H}-= 


k  -  1— >K 


> 


<E 


L(  1  +  k  -»kV- [WS  ♦  TS  -  21^  ♦  10  log  ^mTIoloT^^-^SL^ 


Not*  6 


TN  ♦  10  log  Af  -  H.    ♦  RD  _>M(TN) 


k 


N@t*  7 


arc  sin  (D^/300k)  ^^ 


(&&  i 


-  10  log  bjC^  -  |i)  - 

>  10  log 
179 

b2(0k  -  ji) 
&3(p)  b^(p  - 

-  10  log  (eosfc^  -  p )/«©»&) 
■  200k  o^/c)-^^). 

-  10  log  (1/360)  ^ 
P  = 

1 
■  -  180 

(Note  9 

WS  -  2Hk  ♦  10  log  «fl 

"  ''A 

♦  10  log  (c 

AT/2)  +  10 

log  100k  ♦  RD 
>M(RS)k 

M(SN)  ->  A 


M(TN)k— >B 


RJ2 


D->A 
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M(RS)k-^>B 


RJ2 [D 


D-^M(BG), 


fttote  11 


1 L 


HEL-5> 


^ 


n 


kr 


SL^-M(BG)k^Bjk 


k  -  1-^K, 


OUTPUT 


jkt  J  -  "09  j,  «J9  K  -  19  9    A, 


Note  Lj  This  flow  chart  neglects  the  effects  of  ambient  sea  noise 

of  volume  reverberation.  Both  are  considered  to  have  negli 
gible  effect  on  the  interference  background  masking  level 0 
If  it  is  desired  to  consider  the  volume  reverberation,  it  i 
M(Rv)k  =  WS  -  2^  +  10  log  mv  -(j^  +  10  log  (-£-£!) 


+  20  log  (100k)  +  RD;  k   1,  •••,  K 

where  M(R  )  =  volume  reverberation  masking  level  at  range  100k  yards  (db 
N  v  k 
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WS  =  on-axis  source  level  (db  at  1  yard  from  transducer) 

H,  as  one  way  transmission  loss  at  range  100k  yards  (db) 

10  log  m  =  volume  scattering  coefficient  (db) 

(j  ).  =  volume  reverberation  index  of  the  transducer  (db) 
v'k 

c  =  velocity  of  sound  in  the  medium  (yd/sec) 
AT  =  ping  length  (sec) 

RD  =  recognition  differential  of  transducer  (db) 
The  index  (j  ),  may  be  approximated  by 

v  K 

(jv)k  =  -io  log  ^L  Jj^O)  b2(j)  b3  (p)  bh  (  p°20Qkms)d^. 

iC  ~  1.   000o  & 

where  o)_  =  search  rate  of  turn  (deg/sec) 

3 

j  =  off°axis  angle  in  the  vertical  plane  (deg) 
P  =  off-axis  angle  in  the  azimuth  plane  at  time  of 
pulse  transmission  (deg) 
b..,  b  ,  b  ,  b,  ss  normalized  pattern  functions  of 

the  transmitting  (1,  3)  and  receiving 
(2,  k)   patterns  in  the  vertical  (1,2) 
and  azimuth  (3,  ^)  planes 
_Q_  =  solid  angle  subtended  by  the  sound  energy  beamQ 

This  in  effect  assumes  a  medium  of  infinite  extent 
in  the  direction  of  the  sound  beamc 
Note  2s   M(SN)  =  SN  +  10  log  A  f=DR  +  RD 

where  M(SN)  =  self-noise  masking  level  (db) 

SN  ■  equivalent  isotropic  self-noise  level  at  search  speed 
and  depth  (db) 
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A  f  =  receiving  bandwidth  (cycles/sec) 

DR  =  receiving  directivity  index  (db) 

Note  3:   2H  ^  WS  +  TS~M(SN) 

where  H  =  maximum  possible  on-axis  range  for  acquisition  (yd)„ 

This  neglects  noise  levels  except  the  self-noise  level 

and  neglects  the  weapon  enabling  level. 

TS  =  target  strength  (db) 

Note  ki     H,  =  20  log  100k  +  100  ak 

where  a  =  attenuation  coefficient  (db/yard)  and  it  is  assumed  that 

the  transmission  anomaly  is  zero. 

Note  5s  This  loop  finds  the  upper  limit  (lOOK)  on  the  acquisition 

ranges  to  be  considered „ 

Note  6:  M(TN),  =  TN  +  10  log  A  f-H,  +  RD^  k«l,  ...s  K 

where  M(TN),  =  target  noise  masking  level  (db) 

TN  =  radiated  noise  level  of  the  target  (db  at  1  yard  from 

transducer) 

DWS 
Note  7:  <\  =  arc  sin   — ^ 

where 

ft.  =  the  angle  above  the  horizontal  at  which  an  arc  of  radius 

100k  yards  and  of  origin  at  search  depth  will  inter- 
sect  the  surface  of  the  water 
D   =  weapon  search  depth  (ft) 

Ww 

Note  8;   (j  ),  =  surface  reverberation  index  of  the  transducer  (db) 

S  K. 

u  =  angle  of  attack  (pitch)  of  the  weapon  (deg) 
This  is  an  approximation  described  in  Chapter  k   of  refer- 
ence  k „ 
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Note  9:  M(RS)i.  =  surface  reverberation  masking  level  at  a  range  of 

100k  yards  (db) 
10  log  m  =  surface  reverberation  scattering  coefficient  (db) 
Note  10:  This  subroutine  adds  the  effects  of  two  noise  levels  in 

decibels.  A  rounding  off  process  is  employed  based  on  Figure  ; 
of  reference  3.  A  difference  in  levels  of  more  than  10  db  pro= 
vides  a  sum  effect  equal  to  the  larger  of  the  individual 
levels.  Two  equal  levels  have  a  sum  effect  equal  to  3  db 
higher  than  either  individual  level. 
Note  11:  WEL,  =  weapon  enabling  level  (db).  Values  out  to  a  range 

of  100K1  yards  are  inputs  where  K  is  estimated  such 
that  it  will  be  larger  than  the  K  found  in  the  program, 
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APPENDIX  C2 
SCANNED  VOLUME  RADIALSS  GENERAL  CASE 


INPUT 


INPUT 


PRR ,  «9  «« ,  D, 


TM 


APPENDIX  CI  OUTPUT 


c 


200  oyb  -^  <PX  I (l-^lV^-j  =**  -  Kh  <P%  ♦   (1  -  1)  »1-^Pi 


P 


Note  1 


i   -   1— 5»l| (0-^>J  > t 


D 


1  ♦  k^>k)-i 


j^l — <i->i)JPpi  -  (116^1)  ♦  l)^  ♦  (i  -  i)01_^pi 
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RJ4 


Xl-^>m> <l-=>k) (0->nV 


I 


fNote  9j 


pjmk  :  PJm(k*l)f 


-(  1  ♦  k  -Vk~) — (  1  ♦  n  -»"n>- ■ 


? 


Not*  10 


Pjmk-^Pj(k-n)H £°°*  '  c/2^}^-lc/2PRR8-^Rj(k^ 


100k -=>R 


J(k-n) 


j  :  •* 


<D?M  -  V^ 


J(k=n)   °  -sin(j*|i)c©sp 


'JOmO 


Note  11 


(DTM-DWS 


■srj^TwSpjjj^  ■*■  i<fc^) 


j(k-B)    '  sin(J+i»)co»pj(k_n)" 


"W> 


^»R. 


»in(J*|i)eo8pj^k_n^   '     *   J(k«o 


fNote  12 


pjo«)_m 


k  -  B  -^JK. 


(RTN 


-{  k  -  1  -9>~kl (n  «  2— =>n)- 


{l  +  k-»k) — (T) 


Pjl(k-n)   *  pJ2k 


.(n  -l-»n) (2^m) (T) 


Note  13 


k  -  1-^k) <2-=>m) U  -  2-»b^ (T) 


100 


-(l-j>a) <f  1-^1>— | 


Note 


10  log  \>3(t±)  -  10  log  b4(pt  ♦(!  =  !)  01)-^A 


m  :  1] 


Ai  ♦  i-^i"V 


P(i  -  i)~^fy 


1  *  k-» 


pi  "^  Pjlk 


<j 


i)l 


k :  i 


k  -  l->: 


J  -  1-^J»« 


j  ♦  1^J» 


(RTN 


Note  1:  The  B  .  are  all  possible  transmission  angles  in  plane  zero 
which  meet  the  requirement  that  the  echo  will  be  received 
in  the  interval  I  -  <*»  «  "J  . 

Note  2%     The  overlap  in  plane  zero  is  computed  first  to  provide  the 
angle  B1  .  Then  B1  -  &  and  p'  define  the  limiting  angles 

of  the  Pjmk' 
Note  3;  A.,  .  ■  A  .  (  P  )  for  the  i  th  P  . 
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Note  kt     After  the  computation  of  K,  which  indicates  that  the  maximum 

J 

acquisition  range  in  plane  j  and  all  of  the  8   „  in  plane  i 

jmk 

have  been  computed,  the  subroutine  ends  0 
Note  5?  This  indicates  that  there  is  no  acquisition  probability  in  the 
plane  j  and  changes  the  range  of  planes  to  be  considered  from 

Note  6 1      In  this  case  the  B  .0.    is  limited  by  the  0  B  and  would  have 

been  larger  if  the  overlap  between  pings  had  not  been  removed 
from  the  considered  scanned  area „ 

Note  7;  This  loop  starts  with  the  leftmost  angle  (  B°  ~(l~l)  ^  ) 
to  be  considered  and  computes  each  angle  to  be  considered 
indexed  by  i  and  incremented  by  0=  s    to  a  maximum  angle  B-r^B0  o 

Note  8:  This  point  is  the  commencement  of  the  loop  which  will  compute 

the  B  .  and  R.,  starting  with  j  ■  =  J„ 
jk       JK 

Note  9 i     There  may  be  several  equal  3  #   because  the  scan  azimuth  limit 
were  reduced  by  the  overlap .  The  first  of  the  equal  p    to  be 

considered  will  then  be  the  one  of  maximum  k0 

C        / 
Note  10;  The  ping  interval  range  is        yards/ping „ 


Note  11:  The  range  at  which  the  radial  at  angle  B  .  ,  intersects 
surface  of  the  water  for  j  5*  -  Ji  is  computed  below  „ 
See  Figure  20.  The  distance  0C  is  desired „  The 
is  (DWS/3)  yards. 

length  (CB)  =  Dws  /  3  sin(j^a) 
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Figure  20.   Computation  of  radial  length  for  angle  B  ,  „ 

The  plane  j  for  which  the  scan  outline  is  shown 
is  perpendicular  to  the  vertical  xy  plane „  Angles 
CBO  and  DBO  are  right  angles  0   Length  A$[3l 
search  depth  of  D 


d 


'A- A 
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length  (OC)  =  length  (OB)/c©Sp\*k  =  D^/3  sin(j+ji)  eosp» 

For  j  s-[i  the  development  is  the  same  except  that  (j  +  fli  ) 

is  negative  and  the  considered  depth  is  the  difference  between 

the  maximum  target  depth  D_w  and  the  weapon  search  depth 0 

TM 


Note  12°   If  to  the  left  of  the  radial  of  maximum  range  K, 

J 

and  the  following  angle  will  have  a  larger  k0 
If  to  the  right  of  the  radial  of  maximum  range  K„  then  m  =  2 
and  the  following  angle  will  have  a  lesser  k.,  There  will  be 
one  or  two  radials  of  maximum  range „  That  is 


PjlKj  —  Pj2Kj 


Note   13s      If     P-iir     =     P  -oif  *      t^e  Pr°8ram  skips     p  401f    „     Otherwise   the 


j  J      j  J      j 

program  includes   it„ 
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APPENDIX  C3 
SCANNED  VOLUMES  J  GENERAL  CASE 


0i 


INPUT 


— (  APPENDIX  C2  OUTPUT) (0-^>l; (RJ5)-i 


Note  1 


OUTPUT 


"  <  9    J    =  J*  »  » 


SV 


01 


pjk»  V  j  =  JV>  00%  J"^sl5  0O%  Kj 

P!»  &«  |*«  <%»  Dyse  DTM8   PRRg 


RJ5) (J«-^>J 


-<1  ->k/- 


0  ->  A. 


-> 


RJkRJ(k^l)8ln(Pj(k.l)  '  Pjk)  *  AJ-^Aj| U-UL^k 


{All±l 


J 


9> 


(2/3)ljH(. Dte^fri))***  1  Unih/2)  *  SVn-^SVil 


{TTT^jj 


Note  It      SVQ.  is  not  required  in  the  program  but  it 
to  print  it  here  as  information 0 


be  desired 
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APPENDIX  C4 


SPECIAL  VOLUME  ELEMENTS:  GENERAL  CASE 
INPUT  INPUT 


D"5  a 


O   O    0 


9  D"5  u  =  29 


0    0    0 


D»S 


2, 


*9°  u  = 


rdu,  MBDVdu,  RBDV^  d 

AWdu;  d  =  D+l,  oo° 
BM7duJ  d  =  D+l0    aoe 
ASV^uJd  =  l,    — 
BST^J  d  =  1, 

qSVQuJ  q  =  A9   B: 

hAds  BAId>  tenW  d 
JA0d9  tan0AOd?  d  ■  19 

lM9  hid9  tanW  d  s  2 
BBOd>  W^d-1,    — 

PP         •■»     —     Tfl  °  °  «         TO 

1k*      Ik9   *  '  9 

*DS9  SWD9  STS9  WOSD9  M1V   M»  U 
PS  ^9  JAv  D^o  DTM9  PRRg 


0    0    0 


o  o  © 


U 
U 


29    °°\  U 

Is  o  e 
0  s 


B 
h. 

B. 


©  ©  i 


\  J»S°  k 


9  D» 


ooe 


9      &.£ 


106 


(Note  11 

sy — (j-»j  ) — |o  ->atj 


Note  2 


1 


H.8in(jnO    I  H/.  ^sinO-l^)"]^ |l  ♦  ATJ-^ATj[ 


5F 


3 


1— (J  -l-a»j)- 


&TJ  ->TJ 


ATJ  -  1  -9-ATJ 


^^> 


Ifetg  3  1 


.  H.cos(j*ji)^>BI0 


Hj»ln(j*|i)-H >HLQ 


0-->H 


(J«'*l) 


L-  H 


( j-ATJ)COfl  <*****>  -»Bq0Q 


I 


Note  4 


<OEF} — (OE?>-] 


©—.[Note  51 


H(j*d-r)sln(J+d"r^)^BSd 
HU*d-r)cofl  (j+d"r^)  "^^d 
H(j-ATJ-d*8)8in(J-ATJad>s>Ji)-^SSd 
H(j-ATJ-d+s)COB  (J"ATJ^*s^)-^SCd 

H<J^TM*.*l),to(^T"*"1*)  ^SSd 
H(.1-ATJ-d+a+l)coa  fj^g^^j  -»  sc£ 


(RJ7 


Note  7 


BC. 


■B 


qld 


BId  *  Ld~^BqOd 


* 


Note  8 


"@— I 


L@- 


I— ^^ (1  »  d^>d~> ( 


1  * 


Kd  *  CT(Bql(d»l)   "  BCd)"^>'BqId 
qld  qOd 


J9y 


W10 


L-^jii © 
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H(j-ATJ-d+3)9ln(ja°A™+3^)  -^d 
H(j-ATJ-d+s )cos  0-AT«*8+|i)  -^BCd 
H(j+d.r)8in(J+^-r*l*)^  SSd 
H(j-»d-r)co8  (J*3"1^  -^  scd 

H(jKl-r-l)co8  O+^-r-1^) 


SC! 


Not©  7 


RJ7 


BC— =*B 


qOd 


BqOd  "  Ld 


B 


qld 


~®H 


MTTT^i^n — (TTZ 


■r 


-^7)- 


Kd  -  CT(»d  -  Bq0(M))^B( 


qOd 


BqOi  "  »  ^Bqld 


(lU^) (RJ1Q)— r 

' — 6jlS) © 


Note  15 


7U 


{  1  ♦  d-^d    ) (0-g»u>- 


H(j+d.r)co8(j+d-r^)-^BqId  ■(j«4^)«ta(J«*"**)-^»d 


i 


Bqld  "  U  TES   •  ° 


t-X  i  ♦  u -^»q  ) — (u  :  u 


0-e>E 


^^qld-^0 
V(d-l)  -uTES-^D 


^^  n    d  («+] 


IC 


i  ♦  u_^>tt~) — c§:  j£r 


Q-^qM? 


d(u^l) 


n 
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2H^n^d(u-D  -  ^au-^svodu' 


L(  1  »  n  -^xi~) (u  :  U 


Bqtte 


tan  ( j  *d-r +p )  -5>  C 


»-M*CMt*i 


L 


tan(j+d-r^)-^C 
B  Id  -  (u-l)TES-^D 


C12) (d  ♦  l-^d~) — (2-»u>-* 

4 1  ♦  n  -»«") — (^:  5 


^dC^l)  -  *■&■♦**. 


Q 


d-^D« 


RTN 


Not*  6 


(RJ7)J I  BCd  -  SCd     -   (BSd  -  SSd)  SCJ  -  SCd  /(SSJ  -  SSd)  ->I»dh 


rRJ8] 


111  (d«i)  "  "a-^V* 

BSd-^HLd 


RW 


TN 


(RJ9 


Bql(d-1)  "  BqM"^bId 


RTH 


109 


B  H-   (u-l)TES   I  0 


A 


L-^mj) (rji^- 


IqDVdu  ♦  ^mdu  +  ftl^^qSV^ 


= (RJ12 


A  -=>D 

tan0qH-^C 

0-=^E 


rjh 


B  -  F  -5*-  LqDV 


du 


I*)Vdu  ♦  ^DVdu  i  Rq^^qSV^ 


R5 


BqO(d-D   "    (U~1)TES    :  ° 


(RJ12} 


B-^ttjDV 


du 


M,DVdu  ^  RqDVdu^qSVodu  {— Q 


BqOd  "    (U~1^TES    •  ° 


BqO(d-l)-    fr^Dro^A 


tan^-^C 


D 


B^.-   (u-l)TES  -^D 
qOd 


A-=>E 


Q |  F  -  B  ^>RgDVdu 
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Note  13 


RJljL-L 


A^tan(b/2)  h  ,^B 


RTN 


fNote  U 


RJ14M- 


(1/3)  tan (6/2)  C(D3  -  E3)-2>F 


RTN 


RJ15} 


B^G 


Bql(d-1)  "  ("-DraS-^A 


RJ13 


G  -  B  -^MjDV 


du 


RTN 


(RJ16 


^H^0 


A-5^D 


B qB  -  (u~l)TES-5>E 


RJU) 


B  -  F-S>LqDV 


du 


RTN 


Note  11 


M- 


HL 


^  Ld  -  L(d^l) 


■^HLj 


d-s>D 


-( i » r-^ry 


(HL^  -  ^d'' ^(d-l)  "  ^d 


^(d-l)  "  HLd^^hqd 


RTN 


(RJ18) (_2_-^uJ>  *  |0^»ASVJ (l^»d) 1 


qSV0d  +  ASV  "^ ASV 


-(  1  ♦  d -^»~&~) — fd  :  d 


L  ASV  -^>  qSV, 


Ou 


-(  1  +  u  -^>~u~) (u  :   U 


Note  lj  RJ  6  computes  volumes  either  above  or  below  the  search 

depth. 
Note  2%     The  following  comparison  is  used  to  determine  the  plane 

j  to  be  used  as  a  reference  plane  for  the  computations., 

Ill 


Note  3:  The  following  four  function  boxes  set  up  reference  constants 

for  the  subsequent  computations  <> 

Note  k°.     The  index  d  is  used  to  identify  the  levels <>  The  index  (d=r) 

identifies  the  current  plane  j  above  the  plane  TJ.  The 

index  (d-s)  indentifies  the  current  plane  j  below  the  plane 

T  J-  &  T  J . 

Note  5s   This  comparison  determines  if  the  level  d  is  to  be  determined 

by  a  plane  j>TJ  or  by  a  plane  j  <(TJ=  Atl)„ 

It  is  only  required  for  depths  at  or  above  the  critical  depth^ 

that  is  for  d  -  D. 

Note  6:  This  subroutine  computes  the  distance  L ,  <>   See  Figure  21 . 

d 

As  shown  H/      v  sin  (TJ  +  2  +  u)  > 

H,     x  sin  (TJ-2  +  ji )  and  ATJ  =  0.  The  distance  to  be 

computed  is 

L,  =  length  AB  =  length  FG  ■  length  FE~  length  GE 
d 


h  -  len*th  «  -     £jg|g     (length  HE) 


gth 

Note  7:   If  the  level  length  L,  is  longer  than  the  distance  TES 

traveled  by  the  submarine  during  a  complete  weapon  search 

turn,  the  critical  depth  is  at  a  level  above  L,  but  below 

L/  ,  ,  N.  Then  the  level  defined  by  TES  must  be  found  and 
(d-1) 

designated  L ,  =  Irr  and  the  original  level  dropped  from  con= 

sideration  until  later.   If  L  is  shorter  than  TES  then  it  is 

above  the  critical  depth  and  its  parameters  must  be  computed 

to  find  qSV  .  ,  the  volume  searched  in  the  d  th  level  on  the  u  th 
n  odu' 

search  turn,   (q  =  A,B). 
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TJ+2 


Search  Depth 


Figure  21:  A  sample  volume  portion  for  computation  of  L , .  The  points 
(TJ-  d)  are  the  radial  endpoints  of  the  corresponding 
planes.  The  level  L_  is  under  consideration  and  dig 
(AB  =  L_)  is  to  be  computed.  The  following  are  equ 


BC3  -  SC3 

BS3  '  SS3 

SS3  "  SS3 

sc'  -  sc 


length  FE 

length  BG 

length  HI 

length  HE 
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Note  8:   See  Figure  22.   Subroutines  RJ8,  9,  and  10  compute  the  distances 

and  angles  on  the  right  portion  of  the  figure  which  have  not 

already  been  computed.  Then  subroutine  RJll  computes  all  of 

the  special  volumes  in  level  L, .  After  RJll  the  level  index 

d 

as  well  as  the  index  r  or  s,  is  increased  by  1.   If  index  s 
is  changed  it  means  that  the  plane  which  determined  the  level 
L,  was  above  plane  TJ.  Thus  d  and  s  must  be  increased  together 
so  that  tht   difference  is  constant.   It  is  their  difference 
that  is  used  for  the  comparison  described  by  note  5° 

Note  9:   This  subroutine  computes  the  scanned  volume  in  level  L, 

(where  d^D)  on  each  search  turn  after  the  first.  The  first 
comparison  is  used  to  find  out  if  the  u  th  search  turn  is  the 
last  effective  search  turn  in  level  L , ,   If  it  is  notp  then 
subroutines  RJ12,  15 >   and  16  are  used  to  compute  the  volumes. 
The  cross- section  in  this  case  looks  like  that  of  the  right 
portion  of  Figure  22.   If  the  u  th  search  turn  is  the  last  of 
the  effective  search  turns  in  level  h,s    then  the  normal  cross° 
section  is  cut  at  some  point  as  in  the  left  portion  of  Figure  22 
In  this  case,  the  point  where  the  cross-=section  is  cut  is 
found  by  one  of  the  subsequent  three  comparisons  and  the  cal= 
culation  of  scanned  volume  is  modified  accordingly. 

Note  10:  This  corresponds  to  a  vertical  through  the  transducer  which 

goes  through  a  point  in  level  L  such  as  point  S  in  Figure  22 „ 

Note  11:  This  is  used  to  compute  the  height  of  the  critical  depth 

above  the  search  depth  (length  OR  in  Figure  22)  and  a  correc= 

tion  factor  to  be  used  to  find  the  constants  for  level  L—  „ 

D 
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B 

c 

'(d-1) 

Ld       A 
I*       0 

\g    s 

p 

u 

. 

I 

Search  Depth 


Figure  22:  Vertical  cross-section  through  levels  Let  and  L , . 

The  right  portion  depicts  the  volume  scanned  on  the  first 
search  turn.  The  vertical  OAB  extends  down  to  the  trans- 
ducer . 
Volumes  corresponding  to  depicted  areas  are 


Area 

HID 

DEI  J 

EJK 

CFG  =  EJK 

ABCF  ^  DEIJ 

HDEK 

ABCG 


Volume 
LqDV 


dl 


MqDV 
RqDV 


dl 


dl 


RqDVdu  4   RqDVdl 


MqDVdu  £   MqDVdl 


qSV 


odl 


qSV 


odu 


The  distances  used  to  compute  volumes  in  level  L,  ar« 

d 


length  NQ  =  TES 

length  BD  =  Bql(d-1) 

length  BE  =  Bqo(d_1} 

length  AH  =  B 


length  HI  -  b 


Id 


length  JK  -  length  FG  -  b 


od 


length  AB  =  h 


qd 


qld 


length  AK  =  B 


qod 


length  AR  *  HL 
length  BR  ~  HL 


length  HK  =  L 


(d-1) 


length  DE  -  l, ,  . s 
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Note  12:  This  is  the  volume  computation  for  BqDV,  .  This  computation 

begins  by  turning  the  volume  upside  down  so  that  level  L,  is 

d 

above  level  L,  .   ,\,     Then  the  computation  is  the  same  as  that 
\  d"  L  ) 

for  qSV.,   (corresponding  to  area  ABDH  in  Figure  22)  required 

later. 
Note  13:   See  Figure  23. 
Note  Iki     See  Figure  2*4- . 

length  A  I  =  D  tan0 

area  IJK  =  D2  tan  ( b  /2 ) 

volume  AJKI  =  (l/3)  (area  IJK)  (length  Al) 
=  (1/3)  D3  tan  (&/2)  tan0 
Note  15:  The  volume  computations  (qMV   )  for  levels  below  the  critical 

depth  are  begun  at  this  point. 

Note  16:  The  volume  computations  for  the  last  level  use  B  _,  for  the 

r  qld 

previous  level  and  the  angle  j  +  d  =  r  +  jji  because  there  is 

no  plane  j  on  which  to  base  the  calculations , 

Note  17:  The  calculations  below  search  depth  are  identical  to  those 

above  after  the  H.  below  are  substituted  for  the  H„  above „ 
J  J 

Note  18:  The  final  section  computes  the  scanned  volume  on  the  u  th 
search  turn  separately  for  above  and  below  search  depth. 
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Figure  23:   Volume  computation  of  a  triangular  block.  Length 
EB  =  A.   Length  BC  ■  A  tan  (b/2)o   Length  ED  -  h 
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APPENDIX  C5 
FIRST  TURN  SCANNED  VOLUME:  i  TH  PING 
INPUT  INPUT 


STC»  RI 


r 


APPENDIX  C4  OUTPUT 


(1/3)Dwscot0DS-^RQ|-| 


9 


DWS  STC 
R      +  — _ 

I      3  ain^c  S 


DS     WO 


♦   (i-1)  STg/feRR| 


>-»J>t 


V/i^iyX 


Rjk  cosOnO->Rjk 


Jjk 


=1)6-^  p 


P  :o^ 


p :  i8o 


180  +   p  _=>.  p 


(Note  1 


p  -  180  -=>  p 


180  -  p— =>p 


^Ro2  +  Rjk2-2RoRjkoosP  ^CRJk 


L  arc  sin(R0/1li)sinp  — s>A 


180  -   p  -  A  ^  B 


Rnsinp/sinAcos  (j+n)  — >R 


Jk 


K  1  ♦  k  ^>k    > 


dim 


Hi*l->1   V— 6 


Note  lj   p  is  now  the  angle  between  R  and  the  radial 
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OUTPUT 

iv  svn;  i  =  1,  •••,  m>i 


lQ,  M,   U,  &,  \ 


d;  u  =  29  °°\  u 


♦  wOSD»  W1V  DWS9  DTM9  ?BRa9   P° 
LADVdu,   MADVdu,  RADV^S  d  =  1,    "\  D99  u  =  29 

LBDVdu,   MBDVdu,  RBDV^;  d  =  1,    —  f  DJ  u  =  29    eoe. 

AM7duJ  d  =  D+l,    '*•,  D'J  u  =  1,    °°\   U 

BM\rdu;  d  =  D+l,  •••,  d»;  u ■  =  1,  08\  U 


U 


ASV( 
BSV. 


rMu;  d  =  i,  ••%  d«;  a  =  2,  ooe 
rMu;  d  =  l,  •••,  d«j  u  =  2S  — 

qSV^J  q  =  A,   B;  u  =  2,    "*,   U 

hAd»  BAId'  WAId'  d  =  *■    "°>^ 
BAOd'  tan*AOdJd=lf    — ,  D 
hBd»BBId'  tan0BIdJ  d  -  1,    — ,  D« 
BBOd>  ^W*'1'    '"'S 


U 

i  u 
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APPEND IX  C6 


SCANNED  VOLUMES  SUBSEQUENT  TO  FIRST  TURNS  i  TH  PING 


INPUT 


1    ) (  APPENDIX  C5  OUTPUT") (  1  -^iy^— |  p°    +   (i~3/2)6-=>p 


S. 


D 


p  :  i8o 


< 


180  -  p  -=>p 


180  +   p  -^.p 


□L 


p  »  180  _>  p 


Note  1 


arc  sin     (Rg/R^sinp  — >  A 180  -  p  -  A^>B-,    /T 


H,  sinB/sinA  ->■  CH 


<X^q> 


D  -=>d< 


D» 


{  1  +  u  -g>  u   ) 


D« 


<EEEiD 


0^>ASV[ (l^d> 


L| B  ^  -   (u-l)TES   :  CH 


ASV  +  qSVQdu-^ASV 


> 


CH  -  B  M  *    (u-l)TES  ^>  L 


u 


(CH)Ztan(6/2)(hqd  -  Ltan0qQd) 


(l/3)tan(b/2)tan0qOd  [(CH?-   (B^^  -   (u-l)TESrJ 
(BqO(d-l)  "   (u-DTES)2tan(&/2)hqd 


^>P 


P  +  LqDVdu  +  MlDVdu  +  ASV  -^  ASV 


Bql(d-1)  "    (U"1)TES    :  CHI 


Note  6 


x> 


/^\L[(CH)2  -   (B  I(d  j»  -   (u-l)TES)2Jtan(6/2)h       *  LqDVdu  *  ASV  ^>  ASV 

\j ,  ^ 
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Note  8 


6 


(CH)nan(b/2) 


hqd  '    (Bql(d-l)   -    (U°1)TES  "  CH)tanVd] 


(l/3)tan(6/2)tan0qia  [(CH)3  -   (B  H  -  (u-l)TES)^] 


+  ASV-=>  ASV 


1  +  d-^>  d 


} — Cti!^- (D 


Note  9 


Bqlfi-i)    (U"2^TES    •  CH 


<  fNotelO 


AST  *  qSVQdu-^ASV 


<D 


(l/3)tan(b/2)tan0qia  [(CH)3  -   (B  H  -  (u~2)TES)3] 


qMT,      +  ASV  — *»  ASV 


9 


Bql(d-1)  "   (u-1)^    rCH]^-,  [Note  11 


(CHrtan(&/2)h  d  -  qM^du  +  ASV  ^>  ASV 


Note  12 


(CH)2tan(6/2)jhqd  -   (Bql(d-1)  -   (u-l)TES  -  CH)tan0qId] 

-   (l/3)tan(b/2)tan0qId[(CH)3  -   (BqId  -   (u°l)TES)3]    *ASV->ASV 


I 


(J) — (  1  ♦  d-»d~) (J 


9 


U  ASV-^>qSVlu| (1  +  u^>u~)--H^«^F= 


("Note  13 
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Note  15~] 
3  V  (2-»u>-| 


ASV„      ♦  BSV.    -^SV„ 
iu  iu  iu 


-( i  ->•  u  -=>u"") — (u  :  u 


L<1^1-»1    )— (1 


Note  16 


OUTPUT 


sviu;  1  =  1,  •••,  mo.;  u  =  i,  oo\ 
r.;  i  =  l,  •••,  m-1 


u 


1     M,  U,  6,  Dyg,  DTM,    p«,  PRRg9  Vlt9  v{ 


OSD 
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Note  Is   p  is  now  the  angle  between  R  and  the  bisector  of  the  ith  p 

Note  2%      In  this  case  the  target  volume  cylinder  is  large  enough  such 

that  the  scanned  volumes  after  the  first  search  turn  are  not 

effected. 

Note  3°   In  this  case  the  target  volume  cylinder  is  small  enough  such 

that  some  of  the  scanned  volumes  after  the  first  search  turn 

must  be  modified.  This  is  accomplished  above  and  below  the 

search  depth  and  for  each  search  turn  by  summing  the  modified  or 

unmodified  volumes  over  the  levels  to  obtain  qSV\  ,  First 

T  iu 

each  search  turn  above  search  depth  is  considered,,  then  each 

search  turn  below. 

Note  ki      In  this  case  the  original  qSV  .   is  not  changed , 

Note  5 :   In  this  case  the  original  LqDV ,  and  MqDV   are  not  changed 

but  RqDV  ,   is  modified  and  must  be  recalculated.  The 
du 

original  cross-section  appears  in  Figure  25 .  Using  formulas 

of  subroutines  RJ13  and  RJ14  the  new  volume  P  is  P  -  volume 

FGED  +  volume  ADEC  -  volume  ABGF . 
Note  6:   In  this  case  the  original  RqDV,   is  zeroj,  the  original  LqDV, 

is  unchanged  and  the  modified  MqDV  ,   is  computed  from  the 

formula  for  subroutine  RJ13. 
Note  7:   In  this  case  the  original  MqDV   and  RqDV,  are  zero  and  the 

modified  LqDV,   is  computed  from  the  formulas  for  subroutines 

RJ13  and  Ik   as  the  difference  of  two  volumes. 

Note  8;   In  this  case  the  contribution  of  level  L,  on  the  u  th  search 

a 

turn  is  zero. 
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Figure  25i   Computation  of  modified  volume  RqDV,,  „   Length  FG  is 

dim 

Bq0(d-1)  "  ^u  '  ^   TESo  Area  SSH  represents  the  original 

volume  RqDV,  .   Line  AF  is  part  of  the  vertical 


the  weapon, 
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Note  9:   The  contributions  of  levels  above  the  critical  depth 

now  been  determined . 
Note  10;  The  contribution  of  level  L,  is  not  changed „ 
Note  11:  The  missed  volume  on  the  previous  search  turn  has  been  reduced 

leaving  less  to  be  picked  up  on  the  unchanged  second  turn. 

The  formulas  are  those  of  subroutines  RJ13  and  RJl4„ 
Note  12:  The  missed  volume  in  both  the  current  and  the  previous 

search  has  been  reduced.  The  formula  is  identical  to  that 

discussed  in  Note  7. 
Note  13:  The  scanned  volumes  indexed  by  one  set  of  q^  u,  and  i  are 

now  computed.  The  next  step  is  to  change  u  and  compute  the 

next  set  of  volumes. 
Note  14:   Scanned  volumes  on  the  i  th  ping  of  all  search  turns  above 

search  depth  are  now  computed.  The  next  step  is  to  compute 

the  volumes  below  search  depth. 
Note  15:   The  increments  of  scanned  volume  for  ping  i,  qSV\  ^  are  now 

computed.  The  next  step  is  to  compute  SV.  „ 
Note  16:  Having  found  SV.  ,  the  index  is  changed  to  find 

(i  +  l)u 


1U" 
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APPEND IX  C7 
TARGET  VOLUME!  i  TH  PING 


INPUT 


1  ) (  APPENDIX  C6  OUTPUT  ) (l-^i> 


(JL^JJ  @ 


L|(R02  -  R12)(sin  2(p»    ♦   (1-D&)     -  sin  2(p°   ♦   (1-2)6)   MWLj* 


L(DTM/2)(RQ2  -  Ri2)(sin  2(p»   +  (1-1)6)     -  sin  2(p'    +   (1-2)6)   )->T7± 


OUTPUT 


j  PQ±,  T7jTi  =  l 


9    "••,  M=l 


sviu«  1  =  1,   —,  „t 


M=l°  u  =  1,    00\ 
^  6,  6«,  PRRs9  vlt>  Wpgp 
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Note  L:  The  problem  is  to  compute  areas  as  depicted  horizontally  and 

vertically  lined  in  Figure  16,  Using  Figure  16#   with  origin 

at  point  C  and  x  axis  toward  point  C,  the  equation  of  the 

target  volume  cylinder  circumference  on  ping  i  is 

(x  +  R  f   +  y2  =  R.2 
N     o'    '     1 

Changing  to  polar  coordinates   (r^    p)  with  the   same  origin 

and  measuring  angles  positive  counterclockwise   from  the 

x  axis,    the  equation  becomes 

r2  +  (2R     cosp    )    •   r=(R2  =R2)     ■     o 
o  r  >    l       o' 

Solving  for  r, 

r  =  R,      M  -   (R  A.)2  sin  2p     =R     cosp 
i  o'    i'  o 

To   find   the  area   for  ping  i,    r(r  d   p    )   is   integrated   from 

p'    +  (i  -  2)6     to   p»  +  (i  -   1)&   . 

Thus  _ 

/p'  +  (i-l)b 
area  ■      /  r     dp 

J  p'  +  (i-2)& 

Let  p'   +   (i-l)6  =  B 

p»    +   (i-2)b  =  A 

Then 

B     n  rB 


f     r2  dp  =  f   \fr±2  -  RQ2)  ♦  2  RQ2  cos2p 

-  2  RJt^V1  -   (RAAJ  sinl  cosp  dp 

OX  ox 

J     (R±2  -  RQ2)  dp  =   (R±2  -  Ro2)(B  -  k) 
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f 

^A 


2  R  2  cos2p  dp  =  R 2 
0  r     r         o 


p  +  (1/2)  s±n2p 


Let  x  =  sinp 

Then       dx  =  cosp  dp 


2  R^  ML-  (Rq^)2  sin2p    cos  p  dp 


-B 


2RoRi 


rD  \/i  i  cr^)2  x2  dx 


Let  sinp  =   (RqA^)  x 

Then  cosp  dp  =   (R  A.)  dx 


'  2  RoRi  J       ^1  ~   (\fet>2*    dx  =  "  2  Ri2  J         cos2P  dP 

a  .  R^  [p  *   (l/2)sin2p 


Therefore  J 

-B 
'A 


f      r2  dp  =   (R  2  -  R42)   (sin2B  -  sin2A)  /  2 
Ja  °  i 
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APPENDIX  C8 

ACQUISITION  PROBABILITIES^   I  PINGS 

PPUT .       ,    TNPTIT _         t k 

1  J (  APPENDIX  C7  OUTPUT  ) (   Aw9  W%  W  } (1— »tt)^E 


H  1  •*•  u-^u~) Cu^To' 


^^^H  «55**>  +  «K«*l)&-  *  >/&-».«, 


-(  1    »   U  -=>  U    V 


1    +   S  ^3    ) (1-»U>- 


O-^g 


MU 


^YI^TVl 


L|giu^gs| — Q 


+  S-3>  S 


) — ( i  +  i^il — cT«  5- 


4  1  »  a  -=>"iT^ CjT? 


U 


(W«/Aw)  -   (W8/aw)^R  (i-»r  >-. 


L  w»  ->  w. 


It 


iwlr-WOSD)PRRS/wl|— ^ll 11-^^ 


[^S 


o^PA(Wlr) 


T^T>£f 


KI^a^^a  — <X±I^D 


l-APA^PA(Wlr) 


KjlZZEEiI)- 


Aw  ♦  w 


l(r-l)"^wlr 


Ci> 


OUTPUT 


PA(wlr}'  Wlr»  r  =  1»    "*'  R 
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APPENDIX  D 

PROBABILITY  OF  HIT 

INPUT 


W2jk»  w2jk*  J       ±» 
W2k*       =  ' 


9       U,_  g      k     ~"      J. 


Jk' 


O  O  O  Tf 

9  9      R 


Gl^¥2'^'  G2(Weapon  speed,  Depth),  BW9  D^ 
STM»  SWS'  VT»  '"V  *VX)»  °9  DMIN9  RA9  Qs)  *: 


RH 


RJ19Vt 


DT 


D 


(2/3)DTM  ♦   (1/3)0^  ^HD 


*DT* 


'CT 


360/(%^Te-| 


I>HDMIN  +   (q*l/2)(DTM-DM 


q 


n 


G2(SWS»  Dq)-^WStq 


X  1  +  q^q~) C$f» 


Note  i 


J_ 


3 


-w^^rv^ 


Kl-^k  y-^l-^y^-  G1(wJL)-^S 


'2jk'     ^  "WA 


L    2  P. 


BW/2  +  VT 


sua  ^/(Si.ta   "  stJ 


1-^h 


<A5E*> 


G2(SWA,    (Dyg  ♦  HD)/2)-»w 


2t 


RA¥2t/(SWA°STM}-^W2 


Lrh/(swa  -  s^-^t^ 

^CT  "  *CPA^STM  "   **DT  = 
t       -  t         +    

*  ^PA^WA 

— >  TA 
A 

*W       ^PA                                 SWA-STM 

Kl^q>-<I^T> 


0  -^Aw2n 


G1(SWA»  Dq)-^W 


q  Atq 


(TA  WAtq  *  TS  WStq)PEq  +  Aw2n^AW2n 


^  1  »  g_^g^ ft 


I  Aw'    — 3»-v' 


'2I1- ^"2n 
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°~^PHjk^W2jk»  w2jk^ 


1    -     (l~h)n     -a    P  (v'         .    W»  '      ) 


1 1  *  j-^j  ) — rJ 


K^>^> 


OUTPUT 


PHjk(w2jk>  W2jk>'  W2jk>  W2^  J  "  Ji  °°°»  V  k  =  V".  K 


v^;  k  =  l,  •••,  k 


RJ19) (l-^»q> 


(2q  -  1)/Q2-^P 


{  i  *  q^»"q~}— | 


Note  1*   The  problem  constants  are  now  determined „  Next  the 


probability  P  .,  is  determined  for 


j  =  If       ' '  f    J^j  k  =  1^ 


K. 
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